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20.  ABSTRACT  (con’t) 

geometrical  theory  of  diffraction  (GTD) .  The  technique  extends  the  moment 
method  to  handle  many  problems  that  cannot  be  solved  by  GTD  or  the  moment 
method  alone. 


In  general,  arbitrary  radiators  located  on  or  near  canonical  shapes  or 
combinations  thereof  can  be  solved  using  the  hybrid  technique.  Electro¬ 
magnetic  parameters  which  the  hybrid  technique  can  solve  for  include  the 
near  and  far  fields,  current  distributions,  impedances,  and  scattering  data. 
In  this  paper,  a  transmitting  and  receiving  monopole  antenna,  located  on 
the  surface  of  a  perfectly-conducting  circular  cylinder,  are  analyzed  to 
determine  the  amount  of  radiated  power  which  arrives  at  the  receiving 
antenna,  after  accounting  for  the  reflection  and  curved  surface  diffraction 
effects  of  the  conducting  cylinder.  . . 


The  purpose  of  this  paper  is  to  present  the  technique  and  to  demonstrate 
some  of  its  facility  and  its  accuracy,  as  evidenced  by  the  close  agreement 
of  its  results  to  the  results  of  an  electromagnetic  computer  analysis 
program  currently  used  by  the  United  States  Air  Force  (USAF)  solving  similar 
problems.  Using  both  computer  programs,  antennas  are  modeled  as  one-quarter 
wavelength  monopoles  and  are  positioned  on  conducting  circular  cylinders 
with  radii  ranging  from  one  to  150  wavelengths.  The  power  coupling  between 
the  antennas  is  computed  by  both  computer  codes  and  graphed  as  a  function 
of  the  separation  between  the  antennas.  The  near-perfect  agreement  in  the 
results  of  these  radically  different  analysis  codes  implies  that  their 
solutions,  in  general,  are. correct. 
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CHAPTER  I 


INTRODUCTION 

The  wide-spread  availability  of  very  large,  rapid  digital 
computer  systems  has  lead  to  the  use  of  two  fundamental  methods  of 
modeling  the  electromagnetic  characteristics  of  wire  antennas  on  or 
near  three-dimensional  metallic  surfaces.  One  technique  is  known  as 
the  method  of  moments  whereby  the  interaction  between  electromagnetic 
fields  and  their  induced  currents  on  the  wire  antenna,  as  well  as  on  all 
parts  of  nearby  conducting  bodies,  modeled  as  wire  grids  or  flat  plates, 
is  computed.  Consequently,  this  approach,  which  computes  the  electro¬ 
magnetic  field  from  a  wire  antenna,  is  significantly  limited  by  computer 
storage  and  is  only  practical  when  the  nearby  conducting  bodies  are  on 
the  order  of  a  few  wavelengths  in  size.  An  alternative  approach  is  to 
use  the  geometrical  theory  of  diffraction  (GTD)  which  requires  that 
more  assumptions  about  the  problem  be  made,  such  as  the  radiation 
pattern  generated  by  the  antenna,  but  may  be  applied  to  bodies  that  are 
arbitrarily  large  in  an  electrical  sense.  Much  work  has  been  done  in 
the  recent  past  to  combine  the  method  of  moments  with  the  GTD  in  a 
hybrid  fashion  so  as  to  overcome  the  disadvantages  of  these  techniques 
when  used  separately. 

The  purpose  of  this  thesis  is  two  fold.  First,  the  idea  of  a 
hybrid  technique  of  combining  the  method  of  moments  with  GTD  is 
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entended  to  account  for  the  mutual  coupling  between  two  monopole 
antennas  on  a  large  circular  cylinder  by  means  of  curved  surface  wave 
diffraction.  Second,  the  results  of  this  technique  are  compared  to 
those  of  an  analysis  program  used  by  the  United  States  Air  Force  which 
predicts  electromagnetic  coupling  between  aircraft  antennas.  The  close 
agreement  between  the  results  of  these  two  different  analysis  techniques 
serves  to  validate  both  computer  programs. 

The  basic  hybrid  technique  used  in  this  paper  was  first  described 
in  the  literature  by  Thiele  and  Newhouse  [1],  There,  the  technique  was 
applied  to  antennas  on  and  near  finite  planar  surfaces,  Wedge  diffrac¬ 
tion  theory  was  combined  with  the  method  of  moments  to  account  for  the 
finite  planar  surfaces.  Ekelaan  [  2  ],  building  on  this  previous  work, 
developed  a  hybrid  technique  for  combining  the  moment  method  treatment 
of  wire  antennas  with  the  GTD  to  account  for  reflection  of  electro¬ 
magnetic  energy  from  the  curved  surface  of  am  infinitely  long  cylinder 
as  well  as  diffraction  from  the  ends  of  a  finite  cylinder.  In  the 
present  paper,  curved  surface  diffraction  theory  is  used  to  extend  a 
method  of  moments  thin-wire  analyses  program  to  account  for  the  propa¬ 
gation  of  electromagnetic  energy  around  the  curved  surface  of  an 
infinitely  long  cylinder. 
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CHAPTER  II 


METHOD  OP  MOMENTS  THEORY 

In  this  chapter,  a  fairly  detailed  description  of  the  moment 
method  theory  is  presented  which  serves  as  a  foundation  for  the  hybrid 
technique  presented  in  Chapter  IV.  Since  this  paper  is  a  technique  for 
extending  a  moment  method  analysis  program,  an  extensive  discussion  of 
the  method  of  moments  Is  considered  appropriate. 

The  method  of  moments  is  a  procedure  for  reducing  an  integral 
equation  of  the  form 

(2-1)  f  l(z*)K(z,z')<iz,--Ei 

-'over  structure 

to  a  system  of  simultaneous  linear  algebraic  equations  in  terms  of  the 
unknown  current  i(z').  It  can  be  used  to  determine  the  current  distri¬ 
bution  on  an  antenna,  based  on  the  physical  properies  of  the  antenna, 
its  configuration,  and  even  possible  environmental  influences  such  as 
the  presence  of  a  nearby  conducting  surface.  Thus  the  traditional 
problem  in  antenna  theory  of  deriving  the  form  of  the  current  distri¬ 
bution  on  the  antenna  Is  solved.  Once  the  current  is  known,  radiation 
patterns  and  impedance  can  be  determined  in  a  straightforward  manner. 

The  specific  computer  code  used  for  the  moment  method  portion  of 
the  analysis  is  based  on  the  thin-wire  program  written  by  Richmond  [3]. 
It  requires  that  all  antennas  and  scatterers  consist  of  thin  wires  or  a 
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grid  of  thin  wires.  This  constraint  allows  one  to  make  use  of 
Pocklington' s  integral  equation  to  describe  the  relationship  between  an 
electric  field  incident  on  the  surface  of  a  wire  and  the  resulting 


surface  current  induced  on  the  wire.  Details  of  the  derivation  of  the 
Pocklington  integral  equation  are  given  by  Stutsman  and  Thiele  [4],  The 
general  form  of  the  Pocklington  integral  equation  for  a  wire  of  length 
L  whose  axis  is  parallel  to  the  2  axis  is  given  by 


where 

e_jkr 

G(z,z')  =  free  space  Green’s  function  = 
r-[a2+<a-*')2]l/2 
a  =  wire  radius 


I(z')  ■  filamentary  line  current  on  the  axis  of  the  wire  equi¬ 
valent  to  the  induced  surface  current  integrated  around 


the  circumference  of  the  wire 

2/7^  /  ■'* 

k  =  wave  number  of  the  medium  ■ 

The  parameters  <sj  and  A  axe  the  angular  frequency  and  wavelength 
respectively,  and  jjl  and  S  are  the  permeability  and  permittivity  of 
the  medium  surrounding  the  wire. 

Equation  (2-2)  is  based  on  several  key  assumptions.  First,  the 
wire  is  assumed  to  be  "thin";  the  wire  radius,  a,  must  be  much  less  than 


a  wavelength.  Therefore,  it  is  assumed  that  only  axially  directed 
currents  (z-axis)  are  present.  Secondly,  the  conductivity  of  the  wire 
is  assumed  to  be  infinite,  and  therefore,  all  induced  current  resides 
on  the  surface  of  the  wire.  Since  the  radius  of  the  wire  is  small. 
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this  surface  current,  J(z',9),  is  uniform  around  the  circumference  of 
the  wire  at  any  given  point  along  the  wire.  Therefore,  by  symmetry 
this  surface  current  can  be  replaced  by  a  line  current  which  is  located 
on  the  axis  of  the  wire.  That  is 
f  2*r 


(2-3)  I(z')  -  J(z\e)ad9 


J  0 


It  is  important  to  note  that  eqn.  (2-2)  is  in  the  form  of  eqn. 
(2-1)  so  that  eqn.  (2-2)  may  be  rewritten  as 

fU  2 

(2-4)  /  l( z* )K( z^tz')dz'  -  -  S^(za) 

J  -L/2 

A  physical  interpretation  of  eqn.  (2-4)  is  that  an  electric  field 

incident  on  the  surface  of  the  wire  at  point  z  along  the  wire  axis, 

a 

will  induce  a  surface  current  equivalent  to  a  line  current  distribution 
I(z')  located  on  the  axis  of  the  wire  over  a  length  L  as  shown  in 
Fig.  2-1. 

At  this  point.  It  is  necessary  to  consider  the  boundary  conditions 
at  the  surface  of  the  wire.  From  classical  electromagnetic  theory,  at 
the  surface  of  a  boundary  on  which  there  is  an  incident  field,  one  must 
account  for  a  transmitted  field  and  a  scattered  field  as  shown  in 
Fig.  2-2.  The  relationship  between  these  three  tangential  fields  at  the 
boundary,  *  =•  z  ,  is 
(2-5)  E1  +  3s  »  E1 


However,  a  time  varying  field  cannot  exist  in  a  perfect  conductor  and 
since  the  wire  is  assumed  to  be  a  perfect  conductor,  the  boundary 
condition  on  the  tangential  fields  at  the  wire  surface  becomes 
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INClOENT  FIELD 


i 


w 


(2-6)  31  +  2s  =  0 

Eq  .  (2-5)  can  be  rewritten  as 
(2-7)  Es  =  -  a1 

Finally,  combining  eqn.  (2-6)  with  eqn.  (2-4),  yields 
L/2 


(2-8) 


-L/2 


l(z* )K(z  ,z* )dz'  -  Es(z  ) 
'  m’  m 


A  physical  interpretation  of  eqn.  (2-8)  is  that  a  line  current  i(z') 
located  on  the  axis  of  a  "wire"  of  length  L  produces  a  field  Es(z)  at  an 
observation  point  z  on  the  surface  of  the  "wire"  as  illustrated  in 
Fig.  2-3. 

The  next  important  step  in  the  moment  method  concept  is  to  divide 
the  wire  of  length  L  into  N  segments  of  length  Azq.  This  allows  one 
to  write  eqn.  (2-4)  as 

r  l/2  n  r 

(2-9)  -E^(zm)J  l(z')K(zn*)dz’«  £  In  /  FjK(zm,z’)dz’ 

J  -L/2  n°*l  *^Az 

where  the  current  function  over  the  wire,  l(z')»  is  now  represented  as 
a  series  of  current  magnitudes  at  each  segment  n,  (ln),  multiplied  by  a 
function  which  describes  the  shape  of  the  current  for  each  n  segment, 

Fr  is  referred  to  in  the  literature  as  the  expansion  function.  For 
example,  if  one  chose  to  make  the  current  constant  for  each  segment, 
thereby  forming  a  current  function  I(z')  along  the  wire  consisting  of  a 
series  of  rectangular  pulses,  Fn  could  be  written  as  an  orthogonal 
pulse  function.  Therefore, 


(2-10)  I(z')  -  Y  IF 

i _ _  n  n 

n-1 
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where 

f 1  for  z '  in  z  * 

(2-11)  F  =  /  n 

n  \  0  otherwise 

For  illustrative  purposes,  let  us  choose  the  orthogonal  pulse 
function  as  the  expansion  function.  This  permits  factoring  constants 
out  of  the  integral  and  eqn,  (2-9)  can  be  rewritten  as 

L/2  N 

( I  v»)K(va')ds'  ■ 

-L/2  n-1 


(2-12)  -  E*(z  ) 

z  m 


N 

Lhl,  K<VV)d!s' 

n=4  JAz  * 


n 


Evaluating  the  integral  on  the  right-hand  side  of  eqn.  (2-12), 


(2-13) 


'A, 


X(zm’zn')  dz’  “  G^Vzn’)  A2’ 


Thus,  eqn.  (2-12)  can  be  written  as 

L/2,  N 


f  !  \ 

(2-14)  -Sd(%)  J  £ 

J  -L/2  n-1 


Y_  :n  0  (VY> 

n=l 

Combining  eqn.  (2-14)  with  eqns.  (2-9),  (2-8),  and  (2-7) 
r  L/2  N 

(2-15)  S®(zm)  m  I(z')K(z,z')dz'  -  £  InC(zm,zn')  Az’ 

-L/2  n=*l 

’  IlG<za'zi)  Az'  +  I20(4,,spAs  ♦  ... 

+  ♦  •"  +  V<VV>  - 
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The  physical  interpretation  of  eqn,  (2-15)  is  as  follows  [4],  The  wire 
has  been  divided  into  N  segments  of  length  As'  with  the  current  being 
an  unknown  constant  over  each  segment  (due  to  our  choice  of  the  expan¬ 
sion  function  F).  At  the  center  of  the  segment,  the  sum  of  the 

scattered  fields  froa  all  of  the  N  segments  is  set  equal  to  the  negative 

of  the  incident  field  at  the  point  z  (eqn,  2-7).  The  incident  field  is 

n 

a  known  field  due  to  a  source  located  on  the  wire  (transmitting  case)  or 
a  source  not  located  on  the  wire  (receiving  case). 

By  examining  eqn.  (2-15),  one  may  notice  that  it  resembles 
Kirchhoff's  network  equations 
N 

(2‘16)  EZmIn'V  »-  1.2.3....* 

n=l 

with  which  most  electrical  engineers  axe  more  familiar  (and  comfortable) , 
where 

(2-17)  Zmn  -  C(z  ,z  *)Az 
an  id  n 

V  -  -  E1(  z  ) 

m  z  a 

Since  C(zm,zn')Az  is  based  on  the  geometry  of  the  wire  and  is  given 
for  all  m,  the  values  for  IQ  can  be  determined  by  n  independent 
equations.  In  addition,  the  boundary  condition  of  eqn,  (2-6)  shall  be 
enforced  at  the  center  of  each  of  the  n  segments.  This  is  referred  to 
as  point-matching  which  is  a  special  case  of  more  general  moment 
methods.  Point  matching  results  in  the  following  system  of  equations. 
(2-19)  I1C(z1,z1')Az'  +  I2G(z1,z4,)Az'  +...♦  INC(«1,^(')  -  -B^Zj) 

IlC(z2,z1*)Az'  +  I2G(z2,z2’)Az’  +...+  -  -S2(z2) 

•  •  •  • 

i1c(zn,z1')Az*  +  ♦•••+  VvV)  “  "Vv 
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which  can  be  written  in  compact  matrix  notation  as 


(2-20)  (Ij  iC(2.in)Al'J  -  l-*i(.n)l 


Making  use  of  eqns.  (2-17)  and  (2-18),  eqn.  (2-20)  can  be  rewritten  asi 

(2-21)  [I  ]  [Z  ]  -  V 
n  mrn  m 

Finally,  the  desired  solution  is  obtained  from  eqn.  (2-21)  symbolically 
ast 

(2-22)  »n)  ■  I*...)'1  t  *„  1 

In  practice,  one  does  not  explicitly  generate  [Z  ]  *,  but  instead  may 

n,n 

solve  eqn.  (2-21)  by  using  several  fairly  standard  matrix  algorithms 
which  utilize  properties  of  the  matrix  such  as  its  diagonal 
symmetry. 

Another  important  concept  of  the  method  of  moments  is  weighted 
residuals.  As  the  reader  may  recall,  the  point-matching  method  enforces 
the  boundary  condition  of  eqn,  (2-6)  only  at  the  midpoint  of  any  seg¬ 
ment  on  the  wire.  At  points  other  than  the  midpoint,  z^,  of  a  segment, 
eqn.  ( 2-6)  becomes 

(2-23)  E^+E^-B 

where  R  is  referred  to  as  the  residual.  Thus,  when  R  does  not  equal 
zero,  the  boundary  condition  is  violated;  the  more  severely  this 
boundary  condition  is  violated,  the  greater  will  be  the  chance  that  the 
final  solution  for  the  currents  will  be  incorrect.  By  combining  eqn, 
(2-23)  with  eqns.  (2-8)  and  (2-9) 


(2-23)  R 


L  Xn  +  S~  <%) 

ml  JAz  m 

n 


u 

! 
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and.  for  orthogonal  pulse  expansion  functions, 


N 

a-  T, 

n=l 

In  the  method  of  weighted  residuals,  a  current  function  I(z')  is  chosen 
so  as  to  force  the  residual  to  zero  in  an  average  sense  [4],  This  can 
be  mathematically  expressed  as 


(2-24) 


-L/2 

-1/2 


W  R  dz  -  0,  a  -  1,2,3,. ..,N 
a 


where  W  is  referred  to  as  a  weighting  or  testing  function.  Substituting 
01 

eqn.  (2-23)  Into  eqn.  (2-24)  gives,  in  generali 


and  for  the  orthogonal  rectangular  expansion  function, 


(2-26)/"  (v.  •  ZV<W>M  k*f 

J  -L/2  n-1  J  • 


L/2 


L/2 


W  •  Ei(z  )dz  -  0 
a  z  nr 


The  point-matching  method  of  enforcing  the  boundary  condition  at 

segment  midpoints  amounts  to  using  the  dirac  delta  function 

(2-27)  W  -  *(z  -  z  ) 
m  ra 


where 


and  eqn.  (2-26)  reduces  to  the  expressions  in  eqn.  (2-13). 


In  the  previous  discussion,  the  orthogonal  pulse  and  dirac  delta 


function  were  chosen  for  the  expansion  and  the  weighting  function  so  as 
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to  keep  the  mathematical  expressions  as  simple  as  possible.  However, 
experience  has  shown  that  it  is  desirable  to  choose  an  expansion  function 
which  will  closely  resemble  the  anticipated  form  of  the  current  and  to 
use  that  same  function  as  the  weighting  function.  When  the  expansion 
and  weighting  functions  are  chosen  to  be  the  same,  the  procedure  is 
referred  to  as  Galerkin’ s  method.  Since  the  current  distribution  on 
thin  wire  antennas  is  generally  sinusoidal,  it  follows  that  the  use  of 
sinusoidal  expansion  and  weighting  functions  leads  to  high  accuracy 
with  many  fewer  segments  required  than  if  the  point-matching  method  is 
used.  Fewer  segments  results  in  fewer  segment  currents  (I  )  and 
therefore,  a  fewer  number  of  independent  equations  to  solve.  This  in 
turn  yields  a  significant  savings  in  computer  memory  and  execution  time 
required. 

The  thin-wire  program  written  by  Richmond  is  based  on  Galerkin* s 
method  utilizing  piecewise  sinusoidal  expansion  and  weighting  functions. 
This  is  the  method  and  computer  program  upon  which  this  thesis  is  built 
and  which  extends  the  class  of  problems  which  the  computer  program  can 


handle. 


CHAPTER  III 


GEOMETRIC  THEORY  OF  DIFFRACTION  - 
APPLICATION  FOR  A  SMOOTH  CONVEX  SURFACE 

The  purpose  of  this  chapter  is  to  introduce  a  high  frequency 
method,  the  geometrical  theory  of  diffraction  (GTD) ,  which  is  applicable 
to  bodies  that  are  electrically  large.  This  chapter  will  then  discuss 
the  use  of  GTD  to  determine  the  propagation  of  electromagnetic  energy 
over  the  curved  surface  of  a  circular  cylinder. 

GTD  views  the  propagation  of  electromagnetic  energy  as  rays 
(analogous  to  rays  of  light)  which  are  subject  to  scattering  and 
diffraction  from  specific  parts  of  a  body  such  as  flat  surfaces,  curved 
surfaces,  or  sharp  edges.  Unlike  the  method  of  moments,  GTD  only 
requires  detailed  information  about  the  Interacting  body  at  points  of 
reflection  and  diffraction.  Thus,  unlike  the  method  of  moments,  the 
complexity  and,  magnitude  of  the  analysis  does  not  Increase  as  the 
electrical  size  of  the  interacting  body  increases.  However,  GTD  relies 
on  the  use  of  more  assumptions  than  the  method  of  moments j  primarily, 
the  ray-like  behavior  of  electromagnetic  energy,  which  tends  to  render 
this  method  useless  at  very  low  frequencies. 

The  remainder  of  this  chapter  will  focus  on  the  diffraction 
which  occurs  at  the  surface  of  a  perfectly  conducting  circular  cylinder. 
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In  conventional  GTD  analysis,  the  total  field  exterior  to  a 
curved  surface  will  consist  of  incident,  reflected,  and  diffracted 
rays  which  may  be  divided  into  five  separate  regions  as  shown  in 
Fig.  3-1  [5)«  The  shaded  Region  II  in  the  vicinity  of  the  shadow 
boundary  is  a  transition  region  which  divides  the  lit  zone  from  the 
shadow  zone.  Very  close  to  the  surface,  Region  II  is  subdivided  into 
Regions  IV  and  71  in  the  shadow  and  lit  zone  respectively.  Region  V 
is  a  subdivision  of  Region  III  which  is  in  the  immediate  vicinity  of 
the  surface.  More  specifically,  Regions  IV  and  V  are  close  to  the 
portion  of  the  surface  which  is  a  caustic  of  the  surface  diffracted 
rays 1  whereas,  Region  II  is  in  the  vicinity  of  the  point  of  grazing 
incidence  (Q^)  which  is  a  caustic  of  the  reflected  ray.  Regions  IV, 

V  and  VI  are  therefore  commonly  referred  to  as  the  caustic  or  surface 
boundary  layer  regions.  The  curved  surface  diffraction  hybrid  analysis 
developed  in  this  paper  will  be  exclusively  for  this  surface  boundary 
layer. 

The  theory  and  equations  which  describe  electromagnetic  behavior 
In  the  surface  boundary  layer  of  an  arbitrary,  perfectly-conducting 
convex  surface  are  explained  by  Pathak  and  Wang  [6],  For  the  sake  of 
brevity,  this  chapter  will  simply  present  the  equations  used  to  predict 
the  electric  field  strength  in  the  surface  boundary  layer  from  an 
Infinitesimal  electric  current  (or  current  moment)  which  is  also  located 
within  the  surface  boundary  layer.  According  to  equ.  (31)  in  Ref.  [8], 
the  electric  field  strength  from  a  source  at  point  P*  which  is  observed 
at  point  F  (see  Fig.  3-2)  is  given  by  the  following  expression! 


UJ> 

Q< 

z 


txl 

o 

<r 

Z) 

O  _ 

(O  o 


as 

s!7. 

c  IA 


o 

Nl 

LU 

Z  J 

O  o 

fsj  Q 

t  5 


z  z 

H-  H* 


.  \  \  ' 


i#^ 

£&\i 


o  o 

Q.  CL 

Q  Q 
-1  -1 
LU  LU 

iZ  u. 


or  a. 


uj£  yN 

2<> 

=>  ter 
co  ± 


Fig.  3-1.  The  rays  and  regions  associated  with  scattering 
by  a  smooth  convex  cylinder 


0-1)  2„e(P)  -  S 


where 


({*■„<*  .yi*y2)  -  V*  ■yi'y2)  +  (j£)  Fs(f 

zq  =  impedance  of  the  medium  =*  ^7  ”  //*/£ 

k  =  wavenumber  of  the  medium  =*  2  #"/  A  *  ^  \/ //c& 

s  =»  length  of  the  surface  ray  geodesic  path  from  P*  to  P 

a 

z?  =  radius  curvature  of  surface  along  the  ray  *  - 5 — 

/ S  sin^S 

S  =  angle  of  path  with  respect  to  the  cylinder  axis 
a  =  cylinder  radius 


G(ks}  ^  j  2rs  zQ 


-1 

y^  =  m  kd1 

2  "I  T  J 

y  =  m  kd2 

=  height  of  source  point,  P  \  above  cylinder  surface 

d2  =  height  of  observation  point,  P  .  above  cylinder  surface 

T  =  Torsion  Factor  -  cos  5  for  a  circular  cylinder 
o 

D  =  spatial  factor  «  1  for  a  circular  cylinder 
F  (f  ,y1 ,y2)  =»  "hard"  type  surface  Fock  function 
F  ( £  ,y,y2)  =  "soft”  type  surface  Fock  function 

Pe  =  strength  of  a  current  moment  source  at  P* 
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n'  =*  unit  outward  normal  vector  at  source  point  P* 
n  =  unit  outward,  normal  vector  at  observation  point  ? 

For  £  4  the  Fock  function  may  be  approximated  by  a  Taylor  series 
expansion  as  follows t 

(3-2)  Fh(f  ,y1,y2)  =  V(g  )  -  j  £  _1  r1(  f  )  [  yj  +  y^J 

(3-3)  Fs(f  ,ylfy2)  -  u(  f )  +  |  [u*( f  )  -  \  §  -1  «(f  )]  [y*  +  y|] 

According  to  Pathak  and  Wang  [8],  when  f  is  large,  say  |  > ,6,  a  rapidly 
converging  residue  series  representation  for  the  Fock  functions  may  be 
utilized  as  followsi 

(3.4)  fl/*£  (Tn.)-i  .-J*V 

n»l 

(3-5)  »(  ? )  -  e+J  ^  f  3/2  £ 

n*l 

(3-6)  vl(f)  e*JfV 

n=»l 

(3-7)  u'(f  )  -  eJ,7/4  3/F  f  1/2  £  (1  -  jf  ?rn)e~J?Tn 

n=»l 

where  T  and  T  *  are  zeros  of  the  Airy  function  w0( T  )  and  its 
n  n  2V  n' 

derivative  w2'(T  ),  respectively,  which  are  tabulated  in  Ref.  [8]  and 
repeated  here  for  reference  in  Table  3-1,  Tne  Airy  function,  w2(T),  is 
expressed  asj 

(3-8)  w2  =»  ^7=  j  dz  exp  (Tz  - 

Tz 


w 


TABLE  3-1 


Zeros  of  W ^ (T) 
|r_  I  and  r 


n 


and  W?* 

■  -  lv 


(T) 


-i  /3 


r 


Tn 

V 

2.33811 

1.01879 

4, 08795 

3.24819 

5.52056 

4. 82010 

6.98661 

6.16331 

7. 9^+13 

7.37218 

9.02265 

8.48849 

10.0402 

9.535^+5 

11.0085 

10.5277 

11.9300 

11.4751 

12.8288 

12.3848 

where  integration  contour  Soes  froB1  “°°  to  0  along  the  line  arg  (z)  * 
+27T/3  and  from  0  to  80  along  the  real  axis. 

If,  on  the  other  hand,  £  is  a  small  positive  number  less  than  or 
equal  to  ,6,  one  may  employ  a  small  argument  asymptotic  expansion  for 
the  Fock  functions  as  follows i 

(3-9)  v(  f  )  =  1  -  q.  f  ,  J  fa  3  ♦  .i  »/*  .... 

(3-10)  f  )  -  1  -  f  3/2  +  j  3  +  ^  e-J  9/2.  ... 

(3-11)  ,l(  J)  -  1  ^  €  3/2  -  0  ^  f 3  -  ^  e-J  ’A-  . . . 

(3-12)  ..(f)  3*"V/2  *  J  e-J7r/V/2-... 

In  the  expansion  equations,  eqns.  (3-^)  thru  (3-7),  the  first  ten 
terms  in  the  summation  may  be  used.  For  the  small  agreement  asymptotic 
expressions,  eqns,  (3-9)  thru  (3-12),  the  first  three  terms  may  be  used 
as  shown  in  [8], 

At  first  glance,  the  mathematics  for  computing  the  observed 
electric  field  strength  may  appear  to  be  profuse  and  overwhelming. 
However,  by  confining  the  convex  surface  to  that  of  a  circular  cylinder, 
the  electric  field  strength  of  eqn.  (3-1)  becomes  a  function  of  the 
location  of  the  source  and  observation  points  P*  and  Pj  the  location, 
orientation,  and  radius  of  the  circular  cylinder;  and  the  magnitude  and 
direction  of  the  current  moment  at  point  P',  pg. 


CHAPTER  IV 


THE  HYBRID  TECHNIQUE 


The  hybrid  technique  seeks  to  combine  the  moment  method  with  GTD 
in  an  optimum  manner  so  as  to  yield  an  analysis  procedure  for  antennas 
near  large  electrical  bodies  which  is  both  efficient  and  accurate, 
relying  on  a  minimum  number  of  assumptions. 

The  basic  technique  to  be  described  in  this  paper  was  first  pre¬ 
sented  in  the  literature  by  Thiele  and  Newhouse  [  l].  The  approach  is  to 
first  model  the  wire  antennas  using  the  moment  method  and  then  modify 
the  generalized  impedance  matrix  to  account  for  the  effects  of  nearby 
conducting  bodies  via  GTD.  This  differs  from  other  hybrid  approaches 
which  use  the  moment  method  to  extend  GTD  [7], 

The  basic  moment  method  equation,  eqr..  (2-4), 

f  L/Z 

(4-1)  /  l(z')K(z  ,z’)dz'  -  -  E*(z  ) 

/  EQ  Z  la 

J  -L/2 

predicts  the  electric  field  strength  E^  at  point  ?n  wr.ich  \z  due 

to  the  current  function  I(z')  over  a  wire  of  length  L.  Thru  the  use  of 

expansion  functions,  Fn,  the  wire  can  be  divided  into  N  segments  as 

exDlained  in  Chapter  II  and  so  that  the  field  at  z  can  be  written  as« 

re 


-  Ei(z  ) 
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Applying  a  weighting  function,  Wffl,  to  eqn.  (4-2)  in  order  to  satisfy  the 
boundary  conditions  in  an  average  sense  over  the  length  of  the  wire 
results  in 


r  k/2  N  ^  r  L/2 

W  I  F  K( 2  , z‘ )dz'dz  =  /  w  •  -E1(z  )dz 

ra  n  n  m  /  ra  z'  nr 

-L/2  n=*l  J&z  J  -1/2 


(*-3) 


Interchanging  the  order  of  integration  and  summation 


Following  the  notation  used  by  Thiele  and  Nevhouse,  eqn.  (4-4)  can  be 
expressed  as 


N 

T.  In<V«V>  -<v 

n=l 

where  L(Jn)  represents  a  linear  operator  which  relates  the  surface 
expansion  currents  to  the  field  at  point  m.  Specifically,  it  can  be 
thought  of  as  the  magnitude  of  the  field  at  point  m  due  to  a  unit  test 
current  at  point  n.  For  wires  in  free  space  with  no  other  conducting 
bodies  nearby, 


(4-6) 


L<Jn) 


FnK(Vz,)dz‘ 


ra 

Referring  to  eqn.  (2-16),  eqn.  (4-5) 


can  be  rewritten  as 


(4-7) 


N 

I  Z 
n  mn 

n*»l 


V 

n 


where 


(4-8)  Z  =  <W  ,  L(J  )> 


an 


(4-9)  '/  =  <W  .  -  3^(z  )> 


z  a 


Since  the  linear  operator,  L,  relates  the  expansion  currents  to  their 


electric  fields,  S,  eqn.  (4-8)  could  be  expressed  as 


<4-10)  Z  =  <W  ,  (aS)> 
nvn  a 


where  a  is  a  complex  scalar.  If  conducting  bodies  are  in  the  vicinity 
of  the  wire  antennas,  eqn.  (4-8)  can  be  rewritten  as 


(4-U)  zm-  ,  <Wm.  L'(ja)> 


(4-12)  Zm'  -  <»„,  J?!  4  tS2> 
where 

E^  =  field  arriving  at  observation  point  m  directly  from  Jq 

Eg  =  field  arriving  at  observation  point  indirectly  from  Jn  due 

to  reflection  or  diffraction  from  a  conducting  body 

a  =  complex  scalar  for  E^ 

b  =  complex  scalar  for  Eg 

Using  properties  of  linear  functions, 

(4-13)  Z  '-<W,aB1>+<W,W,> 
mn  nr  1  m'  2 

(4-14)  Z  '  *  Z  +  Zg 
mn  mn  mn 

where 

Z  =  the  "direct  Impedance"  matrix  term  which  relates  the  field 
mn  r 

at  segment  m  that  arrives  directly,  to  the  current  on 
segment  n 
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Z ®  ■  the  “delta  impedance"  matrix  term  which  relates  the  field 

jfln 

at  segment  m  that  arrives  indirectly,  due  to  diffraction  or 

reflection,  to  the  current  on  segment  n 

Z  '  =  the  net  impedance  matrix  term 
mn 

The  hybrid  moment  method  procedure  can  thus  be  summarized  by  the 
following  steps, 

1.  Completely  describe  the  geometry  of  the  problem;  location  and 
size  of  all  wire  segments  and  any  conducting  bodies  in  the  vicinity. 

2.  Determine  the  net  impedance  matrix  terms,  2^'t 

(a)  Determine  the  average  field  at  segment  m,  if  any,  which 

arrives  at  directly  from  a  test  current  of  one  amp  at 

segment  n  (Z  ), 
mn 

(b)  Determine  the  field  at  segment  m,  if  any,  which  arrives 

Indirectly,  due  to  reflection  by  a  conducting  body,  from 

a  test  current  of  one  amp  at  segment  n  (Zr  ). 

mn 

(c)  Determine  the  field  at  segment  a,  if  any,  which  arrives 
indirectly,  due  to  diffraction  by  a  conducting  body,  from 
a  test  current  of  one  amp  at  segment  n  ( 2^n) , 

(d)  For  each  pair  of  wire  segments  n  and  m,  sum  the  direct 
and  delta  impedance  matrix  terms  to  obtain  the  net 
impedance  matrix  term  for  that  wire  segment  pair, 

3.  Define  a  field  source  which  will  generate  an  incident, 
average  field,  V  ,  directly  on  each  segment  m.  (For  transmitting  and 
receiving  wire  antenna  problems, source  field  [V  ]  would  be  non-zero 
only  for  the  transmitting  antenna  segments.) 
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4.  Solve  the  matrix  equation  [  I  ]  *  [Z  ]  *  [  V  ]  to  obtain  the 

n  mn  n 

current  for  every  segment  of  all  the  wire  antennas. 

c.  Use  standard,  simple  electromagnetic  relationships  to  deter¬ 
mine  such  parameters  as  the  antenna  input  impedance,  radiated  power, 
antenna  efficiency,  field  patterns,  mutual  coupling  between  two  wire 
antennas,  etc.,  based  on  the  now  known  current  distribution  on  the  wire 
antennas. 

In  order  to  implement  this  algorithm  in  a  computer  program  which 
would  account  for  cylindrical  surface  diffraction,  the  first  step  was 
to  obtain  an  efficient  moment  method  program.  Richmond's  thin  wire 
moment  method  program  [8],  incorporates  the  Galerkin  method  by  using 
piecewise  sinusoidal  functions  for  both  the  current  expansion  and  the 
weighting  functions.  This  allows  the  program  to  converge  using  roughly 
an  order  of  magnitude  fewer  wire  segments  than  if  a  pulse  basis  function 
point -matching  moment  method  code,  solving  the  same  problem,  is  used. 
Furthermore,  Ekelman  had  augmented  Richmond's  code,  giving  It  the 
capability  to  deal  with  reflections  from  the  curved  surface  of  a  cylinder 
as  well  as  diffraction  from  the  sharp  edges  of  a  truncated  cylinder,  by 
means  of  using  GTD  relationships.  Therefore,  the  program  was  already  a 
hybrid  moment  method  code.  In  both  the  basic  Richmond  code  and 
Ekelman' s  subroutines,  the  Impedance  matrix  terms,  Z^,  and  delta 
Impedance  matrix  terms,  Z^,  are  determined  by  placing  test  currents  on 
segment  n  and  determining  the  resultant  field  at  segment  n.  The  test 
current  distribution  (expansion  function)  is  as  shown  in  Fig.  4-1  with 
n-1  modal  currents  for  each  n  segment  antenna.  Thus,  a  four  segment 
antenna  would  have  three  modal  currents  with  the  end  segments  having 

one  modal  current  each  and  the  other  segments  each  having  two  modal 

2? 


currents.  The  amplitude  of  each  modal  current  is  taken  to  be  one  amp. 

A  typical  segment  n,  which  is  not  at  either  end  of  the  dipole  antenna, 
will  therefore  contain  two  modal  currents  as  shown  in  FI*, 
program  then  determines  the  tangent  field  observed  at  various  points  F 
which  lie  on  segment  m.  This  field  at  P  is  the  vector  sum  of  the  field 
contribution  from  modal  current  I  (ETl)  and  modal  current  II  (ET2). 
Numerical  integration  of  the  field  contributions  is  performed,  resulting 
in  the  impedance  matrix  term  or  delta  impedance  matrix  terra  for  segments 
m  and  n. 

For  computing  the  surface  diffraction  impedance  matrix  terms, 
sd 

Z^,  a  test  is  performed  to  determine  if  the  surface  of  the  cylinder 
obstructs  the  "view"  between  the  test  segment  n  and  the  observation 
segment  m.  If  it  does  not,  such  as  when  segments  m  and  n  are  on  the 
same  monopole  antenna,  the  tangential  fields  are  set  to  zero.  Other¬ 
wise  the  tangent  field  components  generated  by  modal  currents  I  and  II 
on  segment  n  are  calculated  by  subroutines  which  incorporate  the  GTD 
equations  as  discussed  in  Chapter  III.  However,  In  this  case,  the  modal 
currents  are  considered  to  be  located  at  the  midpoint  of  segment  n  where 
the  magnitude  of  each  modal  current,  pg  ,  is  found  as  the  average  value 
of  the  current  distribution  over  the  wire  segment  for  each  modal  current. 
It  may  be  computed  in  the  following  manner.  As  Fig.  (4-3)  shows,  the 
eight  segment  antenna  is  modeled  such  that  nearly  one-quarter  wavelength 
of  the  sinusoidal  antenna  current  is  assumed  to  be  distributed  over  the 
seven  radiating  segments.  To  obtain  a  more  resonant  monopole  (almost 
no  reactance  in  the  antenna  input  impedance) ,  the  total  length  of  these 
seven  segments  was  chosen  to  be  ,2428  wavelengths.  Working  in  units  of 
radians  where  one  wavelength  is  equal  to  277"  radians,  the  length  of  the 


Pig,  4-3,  Current  distribution  on  a  1/4  wavelength  monopole  antenna 
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Pig,  4-4,  Segment  modal  current  distribution 
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monopoie  is  equal  to  1.526  radians.  Therefore,  each  segment  is  one- 
seventh  of  this  length  or  .2179  radians  long  on  which  reside  piece- 
wise  sinusoidal  currents  as  shown  in  Fig,  4-4,  These  segment  current 
distributions  axe  normalised  such  that  their  maximum  value  on  the 
segment  is  one.  In  general,  the  segment  current  distribution  can  be 
described  as 

(4-15)  K?)  -  —bwr  sl°  r 

sin  ( - )  / 

'  n  '  ~ 

where 

L  =  length  of  the  antenna  in  wavelengths 

n  =  number  of  radiating  equal-length  segments  in  the  antenna 
By  using  standard  integration  techniques,  the  average  value  of  this 


current  distribution  may  be  determined  as  follows i 


■<*?*> 


(4-16)  I 


av  z  2  TT  L  \  ,  /  2  7T  L  > 

( - )  sin  ( - )  - 

'  n  '  '  n  '  0 


sin  ?d  j 


Performing  the  integration  yields i 
1 


(4-1?)  I 


/2  7TL, 
cos  ( — — ) 


av  (2 Eli)  sln  (121) 

'  n  '  '  n  ' 


Thus,  for  n  *  ?  and  L  *  .2428,  I  -  ,502  amps.  This,  in  turn,  is  the 
value  that  was  used  for  the  magnitude  of  the  current  moment,  pg  ,  for 
the  surface  diffraction  equations.  The  source  and  observation  segment 
Interaction  for  surface  diffraction  is  illustrated  in  Fig.  4-5.  The 
program  assumes  that  the  currents  on  segment  n  will  launch  two  surface 
field  waves?  one  which  will  travel  clockwise  around  the  cylinder,  the 
other  counterclockwise.  Both  fields  will  arrive  normal  to  the  surface 


mm* 


of  the  cylinder  at  segment  m,  be  added  together  as  vectors,  and  dotted 

with  the  direction  of  the  mth  segment  on  the  receiving  antenna  to  yield 

the  total  tangential  field.  As  eqn.  (3-1)  requires,  the  program  also 

dots  the  a  direction  of  the  source  segment  modal  currents  on  segment  n 

with  the  outward  normal  unit  vector  of  the  cylinder.  For  monopoles 

normal  to  the  cylinder  surface,  these  dot  products  are  all  unity.  The 

tangential  fields,  ETl  and  ET2,  are  then  integrated  in  the  same  manner 

as  for  the  subroutines  which  determine  the  impedance  matrix  terms,  , 

£ 

or  reflection  delta  impedance  matrix  terms,  Z _ , 

an 


X- 
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CHAPTER  V 

COMPARISON  OF  RESULTS 
WITH  AN  ALTERNATIVE  ANALYSIS  PROGRAM 

This  chapter  will  discuss  an  alternative  electromagnetic  analysis 
program  which  is  promoted  by  the  United  States  Air  Force  (USAF)  and  will 
show  how  well  the  results  of  the  LEAF  analysis  program  agrees  with  the 
hybrid  method  of  moments  (MOM)  program  for  a  variety  of  cylinder  radii 
and  relative  antenna  locations. 

The  Intrasystem  Electromagnetic  Compatibility  Analysis  Program 
(lEMCAP)  is  a  USA  Standard  FORTRAN  program  for  computer-aided  imple¬ 
mentation  of  electromagnetic  compatibility  (EMC)  at  all  stages  of  an 
Air  Force  system's  life  cycle,  applicable  to  aircraft,  space/missile, 
and  ground-based  systems  [9]. 

The  antenna-to-antenna  coupling  model,  one  of  the  many  coupling 
modes  between  equipment  subsystems  which  E2MCAP  can  analyze,  is 
basically  a  geometric  optics  (GO)  program  which  determines  the  severity 
of  interference  that  a  transmitting  antenna  may  unintentionally  produce 
in  other  antenna  receiving  systems  on  the  same  aircraft  due  to  such 
phenomena  as  the  generation  of  higher  order  harmonics.  To  use  the 
antenna-to-antenna  portion  of  IEMCAP,  one  begins  by  modeling  the 
aircraft  as  a  conducting  cylinder  with  a  cone  attached  at  one  end. 

Wings,  modeled  as  Infinitely- thin  flat  plates,  can  be  attached  to  the 
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sides  of  the  cylinder.  Next,  antennas  are  described  by  their  gain 
patterns,  supplied  by  the  User.  Up  to  three  different  quantized 
levels  of  antenna  gain  may  be  specified  for  the  spherical  sectors  wnich 
enclose  the  antenna.  The  antennas  are  then  positioned  and  oriented  at 
their  respective  locations  on  the  cylinder  that  represents  the  aircraft 
fuselage.  The  power,  frequency  range,  signal  modulation  character¬ 
istics,  and  relative  harmonic  levels  of  transmitters  are  enumerated  by 
the  User  as  well  as  the  antennas  to  which  the  transmitters  are  connected. 
Similarly,  the  power  sensitivity  threshold,  frequency  range,  and  out- 
of-oand  rejection  of  the  receivers,  as  well  as  their  antennas,  are 
specified  by  the  User.  IEMCAP  will  then  determine  the  magnitude  of 
power  delivered  to  a  specific  receiving  antenna  from  one  of  the  trans¬ 
mitting  antennas,  based  on  the  relative  position  of  the  antennas  on  the 
modeled  aircraft,  the  appropriate  gains  of  the  antennas,  and  the 
frequency  of  transmission.  The  magnitude  of  the  received  power  is  then 
compared  to  the  receiver's  power  sensitivity  for  that  frequency  and 
if  it  exceeds  that  sensitivity,  a  potential  electromagnetic  inference 
(EMI)  problem  is  predicted  by  IEMCAP. 

In  determining  the  magnitude  of  the  power  coupled  to  the  receiving 
antenna,  IEMCAP  takes  into  consideration  the  distance  between  the 
antennas  (free-space  loss),  diffraction  around  the  fuselage  (fuselage 
shading  factor),  and  any  diffraction  off  the  edge  of  the  wing  which  may 
lie  in  the  direct  path  between  the  two  antennas  (wing  diffraction 
factor).  These  three  factors,  all  in  units  of  decibels,  are  added 
together  algebraically  with  the  antenna  gains  (also  in  decibels)  to 
arrive  at  the  power  coupling  factor  between  the  two  antennas  for  a  given 


frequency. 


3< 


The  fuselage  shading  factor  used  by  IEMCAP  is  "cased  on  work  by 
Hasserjian  and  Ishimaru  ( 10  1 .  In  their  analysis,  a  function  is  deri'/ec. 
which  relates  the  propagation  around  an  infinite  conducting  cylinder  to 
that  over  a  flat  plane.  This  function  is  approximated  by  IEMCAP  as 
follows  [ 1 1 j i 

(5-l)  SF=-?7rTT7 

where 


(5-2)  *  -  />,  ej 

C 


*  = 


5.478  x  10_3 
3.340  x  10"3 


for  A  <  26 
for  A  2s  26 


?  = 


'.  5083  for  A  <  26 
. 5621  for  A  >  26 


and 


5F^  -  fuselage  (cylindrical)  shading  factor  (dB) 
a 2  =  radius  of  cylinder  (meters) 

9s  =  angle  around  cylinder  of  propagation  path  (radians) 

^  =  wavelength  (meters) 

Dc  =  distance  of  the  cylindrical  segment  of  propagation  path 
It  is  possible  to  compare  the  performance  of  the  hybrid  moment 
method  against  that  of  IEMCAP  in  computing  cylindrical  diffraction  loss. 
As  part  of  the  standard  output,  IEMCAP  will  print  the  transfer  loss 
between  two  antennas  in  units  of  decibels  (dB),  This  corresponding 
number  can  be  obtained  from  the  hybrid  moment  method  program  in  the 
following  manner.  The  input  power  to  the  transmitting  antenna  can  be 
determined  as 
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(5-3)  Pln  -  Seal  Ik*J 

where 


P,n  =  transmitting  antenna  input  power  (watts) 
y  =  complex  voltage  at  the  antenna  terminals  from  a  voltage 
generator  (volts) 

-ft- 

E  =  complex  conjugate  of  the  complex  current  at  the  antenna 

& 

terminals  (amps) 

The  power  which  is  received,  at  the  other  antenna  is  found  by 
observing  the  complex  current  induced  in  a  matched  load  impedance  of 
36.2  +  j.0104  ohms  (input  or  terminal  impedance  of  a  "near"  resonant 
one-quarter  wavelength  monopole)  connected  across  its  terminals  and 
applying  the  equation 

<5-*) 


where 

*out  =  Power  developed  at  the  receiving  antenna  terminals  (watts) 
Er  =  complex  current  at  the  antenna  terminals  (volts) 

3Tr*  =  complex  conjugate  current  at  the  antenna  terminals  (amps) 

=  impedance  connected  across  the  antenna  terminals  (ohms)  = 
36.2  +  j.0104 

Finally,  the  transfer  loss  or  power  coupling  factor  may  be  computed  as 

(5-5)  FCF  -  10.  log10  |  | 

in 

where 

PGF  =  power  coupling  factor  (decibels) 

PQ  =  power  at  receiving  antenna  (watts) 

P^n  =  power  input  to  transmitting  antenna 
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The  power  coupling  factor  includes  the  gain  of  both  the  transmitting 
and  receiving  antennas.  In  order  to  eliminate  deviation  between  the 
results  of  the  programs  due  to  different  values  of  antenna  gains,  the 
effective  maximum  broadside  gain  of  the  monopole  modeled  by  the  MOM 
program  was  determined  in  the  following  manner,  A  test  case  was  per¬ 
formed  in  which  the  transmitting  and  receiving  antennas  were  located 
parallel  to  each  other  on  the  surface  of  a  ten  wavelength  radius  cylinder 
and  separated  by  a  distance  of  exactly  ten  wavelengths  as  shown  in 
Fig.  5-1.  These  wire  monopole  antennas  had  a  wire  radius  of  ,00005 
wavelengths  and  were  exactly. 2428  wavelengths  long.  The  total  power 
coupling  between  the  transmitting  antenna  and  the  receiving  antenna  was 
determined  by  using  the  procedure  previously  described  and  found  to  be 
-32.65  dB.  Next,  the  far-field  free  space  loss  was  computed  using  the 
familiar  expressioni 


(5-6)  FSL  -  10.  log 


16 


where 

FSL  ■  free-space  loss  (dB) 

D  =  distance  between  antennas  (wavelengths) 

Using  a  value  of  10  for  D,  the  value  of  FSL  is  approximately  -41.98  dB. 
The  difference  between  these  two  values  is  accounted  for  by  the  effective 
gain  of  each  of  the  two  antennas  (assumed  to  be  equal).  Therefore,  the 
antenna  gain  is  found  as 

(<-7>  Cain  ■ 


or 


(5-3)  rjain  -  4.665  dB 


10 \  RADIUS 


Fig.  5-1.  Determination  of  monopole  antenna  gain 


As  a  side  point,  this  is  very  close  to  the  ideal  maximum  broadside  gain 
of  5.16  dB  for  a  resonant  one-quarter  wavelength  monopole  antenna  on  an 
infinite  flat  surface.  The  gain  of  4,665  dB  derived  by  the  MOM  code  for 
the  antenna  on  a  curved  surface  was  used  as  input  data  for  the  TEMCAP 
code  for  the  one-quarter  wavelength  monopole  antennas. 

The  next  step  was  to  choose  an  antenna  configuration  for  both 
programs.  For  the  first  configuration,  a  transmitting  monopole  antenna 
was  fixed  on  the  surface  of  a  circular  cylinder  while  the  receiving 
monopole  antenna  was  positioned  in  the  plane  of  the  transmitting  antenna 
perpendicular  to  the  cylinder  axis,  at  various  angular  separations  from 
the  transmitting  antenna  as  shown  in  Fig.  5-2.  The  radius  of  the 
cylinder  was  initially  chosen  to  be  ten  wavelengths.  This  corresponds 
to  the  fuselage  of  a  cargo  aircraft  such  as  a  C-141  or  a  C-5A  at  around 
one  gigahertz  (GHz),  close  to  the  frequencies  used  by  airborne  military 
navigation  equipment. 

Fig.  5-3  illustrates  the  excellent  agreement  between  these  two 
radically  different  approaches  for  calculating  the  power  coupling  factor 
between  two  antennas  due  to  cylindrical  surface  diffraction.  Here, 
values  were  computed  at  ten  degree  intervals.  The  only  significant 
difference  in  results  between  the  two  methods  is  when  the  two  antennas 
are  separated  by  large  angles.  At  these  angular  separations,  both  the 
short  and  long  path  diffracted  field  contributions  are  similar  in 
magnitude  at  the  receiving  antenna.  However  the  difference  in  the 
phase  angle  of  the  arriving  fields  results  in  a  reinforcement  or 
cancellation  effect.  3y  plotting  only  the  values  between  160  and  130 
degrees,  this  phenomenon  can  be  more  clearly  seen  as  Fig.  5-4  shows. 


POWER  COUPLING  FACTOR  (OB) 

l GO . 00  -120.00  -80.00  -40.00  0.00 


RADTUS-  (i)  =  10.0 
FEMCAP  CODE 

— hybrid  mom  code 


'G.CC  40.00  80.00  120.00  ISO. 00 

ANGULAR  SEPARATION  (DEGREES) 


200.00 
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The  MOM  code  data  in  Fig.  5-^  is  very  similar  to  a  standing  wave  pattern 
on  a  lossy  transmission  line  in  the  vicinity  of  a  discontinuity.  The 
distance  between  null  points  is  one-half  wavelength.  One  observes  that 
very  close  to  180  degrees,  the  total  power  received  by  the  second 
antenna  can  vary  by  almost  as  much  as  20  decibels  for  a  displacement  of 
one-quarter  wavelength! 

Next,  the  radius  of  the  cylinder  was  varied.  As  an  upper  limit, 
a  cylindrical  radius  of  150  wavelengths  was  selected.  This  corresponds 
to  the  fuselage  of  a  C-5A  cargo  aircraft  (30  foot  diameter)  at  a 
frequency  of  approximately  10  GHz  which  is  near  the  operating  frequen¬ 
cies  of  airborne  radar,  precision  approach  systems,  and  electronic 
warfare  (EW)  equipment.  Fig.  5“5  shows  very  close  agreement  between 
the  two  approaches  for  all  angular  separations  less  than  130  degrees. 

As  in  Fig,  5-3,  the  MOM  data  oscillates  for  angular  separations  larger 
than  160  degrees.  However,  since  the  separation  between  the  nulls  is 
less  than  .2  degrees,  these  oscillations  appear  as  a  shaded  area. 

Figs.  5-6,  5-7,  and  5-8  show  that  the  two  programs  generally  yield 
the  same  answers  for  smaller  cylinders  with  a  5,  3,  or  1  wavelength 
radius.  Again,  there  is  some  deviation  for  very  large  angular 
separations  due  to  the  constructive/destructive  reinforcement  discussed 
previously.  There  is  also  somewhat  of  a  deviation  for  very  small 
angular  separations.  This  is  most  likely  due  to  the  geometry  situation 
of  the  antennas  and  cylinder  used  by  the  hybrid  MOM  code.  For  these 
small  angles,  the  electric  field  is  no  longer  propagating  through  the 
surface  boundary  layer  region  described  in  Chapter  III.  Rather,  source 
points  on  the  transmitting  antenna,  as  well  as  the  observation  points 
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Fig.  5-7.  Power  coupling  factor  between  two  l/4  wavelength 
monopoles  oil  a  cylinder  of  radius  3  wavelengtns 
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Fig.  5-8.  Power  coupling  factor  between  two  l/4  wavelength 
monopoles  on  a  cylinder  of  radius  1  wavelength 


on  the  receiving  antenna,  may  more  properly  be  considered  as  points 
well  off  the  cylinder  surface,  especially  for  the  antenna  segments 
farthest  away  from  the  cylinder  surface.  Therefore,  a  different  set  of 
surface  diffraction  equations  should  be  used  which  permits  the  source 
and/or  observation  point  to  be  outside  the  surface  boundary  layer  [6], 

A  case  which  lends  additional  support  to  this  deduction  is  that  of 
the  five-eigths  wavelength  monopole.  The  length  of  the  one-quarter 
wavelength  monopole  was  extended  to  .625  wavelengths.  The  matched  load 
impedance  was  now  found  to  be  73.5  +  j*K)0,7  ohms.  The  antenna  effective 
gains  were  computed  by  the  MOM  code  as  before  and  found  to  be  8.385  dB, 
nearly  equal  to  the  ideal  maximum  broadside  gain  of  8.7  dB.  This  value 
of  8.385  dB  was  used  as  the  gain  of  the  antennas  in  the  IEMCAP  code. 

Fig.  5-9  illustrates  the  results  of  the  two  programs.  Here  the 
deviation  between  the  two  codes  is  becoming  significant.  VJith  the 
antennas  more  than  one-quarter  wavelength  above  the  cylinder  mrface,  many 
more  source  and  reception  points  on  the  antennas  are  no  longer  in  the 
surface  boundary  layer  region.  Thus,  as  for  one-quarter  wavelength 
raonopoles  on  smaller  radii  cylinders,  the  MOM  code  appears  to  be  less 
accurate. 

A  second  type  of  orientation  of  the  antennas  on  the  fuselage  was 
used  to  introduce  torsion  in  the  paths  between  the  antennas.  The 
antennas  now  are  no  longer  in  the  same  plane  perpendicular  to  tne 
cylinder  axis  as  before.  They  are  separated  by  both  an  angular 
displacement  as  well  as  a  displacement  along  the  cylinder  axis  wnicn 
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Fig.  5-9.  Power  coupling  factor  between  two  5/8  wavelength 
monopoles  on  a  cylinder  of  radius  10  wavelengths 


Soecif ically,  two  cases  of  constant  angular  displacements  of  90 
and  180  degrees  were  chosen  where  the  power  coupling  factor  was  computed 
and  plotted  as  a  function  of  linear  separation  parallel  to  the  cylinder 


axis.  These  orientations  are  illustrated  in  Fig.  $-10. 

The  results  of  the  power  coupling  factors  computed  by  both  the 
MOM  code  and  the  IEMCAP  code  for  angular  separations  of  90  degrees  is 
graphed  in  Fig.  5“H.  The  agreement  between  these  two  methods  is 
amazingly  close  for  all  four  linear  separations  computed. 

Fig.  5-12  compares  the  results  of  the  IEMCAP  to  the  MOM  code  when 
the  angular  separation  is  180  degrees.  Although  the  two  codes  produce 
answers  that  differ  by  approximately  6  dB,  this  difference  is  most  likely 
due  to  phase  information  that  the  MOM  code  accounts  for  which  the 
IEMCAP  code  does  not  (as  explained  in  detail  previously  in  this  chapter). 
However,  both  codes  predict  nearly  equal  changes  in  the  power  coupling 
factor  for  corresponding  changes  in  linear  displacement  along  the 
cylinder  axis.  Furthermore,  for  the  180  degree  angular  displacement 
case,  both  codes  show  that  the  power  coupling  factor  will  actually 
Increase  as  linear  displacement  between  the  two  antennas  is  increased. 

Based  on  these  comparisons  between  the  MOM  code  and  the  IEMCAP 

code  in  their  ability  to  determine  the  loss  in  transmitted  power  due  to 

cylindrical  surface  diffraction,  several  conclusions  may  be  drawn.  For 

intermediate  angular  separations  (60  to  120  degrees),  the  deviation  in 

results  between  the  two  codes  is  less  than  one  decibel.  This  is  a  strong 

Indication  that  both  results  are  correct.  For  larger  angular  separations 

(120°  to  180°)  deviation  between  the  codes  can  be  explained  by  the 

reinforcement/cancellation  effects  due  to  phase  information  which  is 

accounted  for  by  the  MOM  code  but  not  by  IEMCAP.  Yet  the  IEMCAP  code 
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values  do  agree  with  the  average  values  of  the  standing-wave  type 
patterns  of  the  MOM  code  results.  Thus,  the  MOM  code  results  appear  to 
be  correct  for  cylindrical  surface  diffraction  involving  large  angles. 

In  addition,  there  is  also  good  agreement  between  the  two  codes  for 
cases  of  linear  displacement  along  the  axis  of  the  cylinder  where 
torsion  is  introduced  into  the  path  between  the  transmitting  and 
receiving  antennas. 

However,  deviation  between  results  of  the  two  codes  is  more  severe 
for  small  angular  separations  (less  than  60  degrees)  between  one-quarter 
wavelength  monopoles  on  small  cylinders  with  radii  on  the  order  of  one 
wavelength.  There  is  also  a  disturbing  amount  of  disparity  between  the 
results  of  the  MOM  and  IEMCAP  codes  when  a  longer  five-eights  wave¬ 
length  monopole  is  used.  It  appears  that  the  MOM  code  fails  (gracefully) 
when  source  points  are  much  more  than  one-quarter  wavelength  above  the 
cylinder  surface  or  when  very  little  of  the  propagation  path  between 
the  antennas  resides  in  the  surface  boundary  region,  as  for  small 
angular  separations  on  small  radii  cylinders.  To  reduce  the  effect  of 
these  problems,  equations  which  account  for  curved  surface  diffraction 
from  source  and  observation  points  outside  the  surface  boundary  layer, 
such  as  described  by  Pathak,  Burnside,  and  Marhefka  [6],  should  be 
considered. 
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APPENDIX  A 


DESCRIPTION  OF  SUBROUTINES 

The  purpose  of  this  appendix  is  to  furnish  a  brief  description  of 
the  subroutines  which  were  added  to  Richmond's  thin-vire  program  in 
order  to  account  for  surface  wave  diffraction  around  a  circular  cylinder. 
For  a  description  of  the  subroutines  used  in  the  original  thin-wire 
program  or  the  subroutines  developed  by  Skelman  which  account  for 
reflection  from  the  surface  of  the  cylinder,  the  reader  is  advised  to 
consult  References  [3]  or  [2]  respectively. 

Subroutine  SGANTC,  listed  in  Fig.  A-l,  is  called  in  place  of 
SCANT.  Subroutine  SGANTC  is  the  same  as  subroutine  SCANT,  described  in 
Reference  [3],  with  two  exceptions.  First  of  all,  subroutine  GGSCYL  is 
called  Instead  of  GGS,  Secondly,  the  statements  which  modify  the  mutual 
Impedance  matrix  C  to  account  for  lumped  impedances  contained  in  the 
ZLD  matrix,  have  been  modified  to  accommodate  the  C  matrix,  stored  as  a 
two-dimension  array  Instead  of  the  original  single-dimension  array  used 
in  Richmond's  program. 

Subroutine  GGSCYL,  listed  in  Fig.  A-2,  is  similar  to  subroutine 
GGS,  of  the  original  thin-wire  program,  described  in  Reference  [  3], 
However,  GGSCYL  will  test  If  the  surface  of  the  cylinder  obstructs  the 
line  of  sight  between  the  end  points  of  the  source  and  observation 
segments.  This  is  done  by  calling  subroutine  VISIBLE  whose  parameters 
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SUBROUTINE  SGAN I C (  IC J . I A . I B , I NM . I  NT , l SC . I  1 . I  2 . 1  3 . JA . JB . MD . N . NO . NM . 

2  NP  .AM.BM,C.CG0.CMM.D.EP2.EP3.ETA,FHZ,GAM,SGD. 

3  X .  1  .  Z.71 0 ,ZS. I2NM.RA0M  ) 

*  *  ,<  ANIC  IS  THE  SAME  AS  SGANT  EXCEPT  THAT  GGSCYL  IS  CALLED  INSTEAD 
C  *  *  01 

COMPLEX  ZG,ZH,ZS.EGD,GD,CGDS,SGDS,SGDT,B01 

COMPLEX  P11.P1Z.P21.P22.Q11. 012. 02  I  . Q2 Z . E P 2 . • P . E TA . GAM . E P 3 

COMP LEX  EPSILA. CWEA . BE  TA , ZARG 

COM'  LtX  P<2, 21.0(2,2)  . CGD(  I NM  ) . SGO ( T NM  ) , C< I C J . IC J  )  . 7L D<  I 2NM  ) 
DIMENSION  X<NP),Y(NP>.Z(NP),D(  1 NM  > . I A(  INM),1B(  I NM  )  , MO i  INM.4) 
DIMENSION  I 1< 1CJ  ) ,  I2l  ICJ  >. I3<  ICJ  )  ,JA(  ICO  )  .  JB(  ICJ  )  ,ND(  I NM  > . I  SCI  I NM  ) 
DAT,,  E0.TP.U0/3.854E-  1 2 , o . 2 8 3 1 8 , 1 . 2566E - 6 / 


2  OS.  MAT  (  3X  ,  ‘AM  *  *  .  E  i  j»  .  3 . 3X  ,  ‘DMAX  »  »  .  E  1 0 . 3 . 3  <  .  *  DM  1  N  --  *.F10.3> 
E  P  -  C  ?  3 

ltC=<N‘N+N>/2 
DO  1  3=1,  N 
DO  10  1  ■=  1  .  N 

10  C  i , 0 ) = i 0 . , 0 . > 

11  .'ONTINUE 

2  S  “  !  .  0  .  .  0  ) 

IP ( CMM.LE .0.  )G0  TO  12 
OMEGA=TP*F HZ 

E  P S I L A  =  CMP  L  X ( E0.  -CMM*1 .Eo/OMEGA) 

CWI  A=i  .0. 1 .  )  ‘OMEGA*  E  P  S  !  LA 

BE  fA=OME GA'SQRT t U0  )*CSQRT< EPSILA-EP  > 

ZARG -BETA*AM 

CAI  L  CBE3<  ZARG . B0I  > 

ZS=8£TA‘B01 /CWEA 
12  ZH=ZS.'  (  TP  *AM*GAM  ) 

DM  I N  =  1 . £  30 
UMAX*-  .0 
DO  20  3= I  ,  NM 
r  =  ;a(3  i 
L  -  I  8  (  J  > 

D'.  3  )  =SQRT  (  (X<K)-X<L>>“2  +  (Y(K)-Y(L>)“2  +  (Z(K>-Z<L))“2> 

IF( 0(0 J.LT.DMIN  )DM1N  =  0<3  ) 

I F ' D( 3  )  . GT . DMAX  IDMAX  =  0<  3  ) 

EGD*CEXP( GAM*D ( 3  )  > 

CGOi 3  >  =  ( EGD+1 .  /EGO  >/2. 

20  SGD  J  )  =  ( EGD-1 . /EGD  1/2 . 

IFCDMIN.LT. 2. ‘AMIGO  TO  25 
IF  v  C\BS< GAM ‘AM ).GT.0.06)GO  TO  25 
Ir(CABS<CAM»DMAX).GT.3.  )G0  TO  25 
1KAM.GT.0.  >G0  TO  30 
25  N  =0 

WRITE ( 6 , 2 ) AM, DMAX ,DMI N 

Rt  f'.'P.N 

30  DO  200  K=1 ,NM 
NO K  =  ID (  K  ) 

PA= I A  < K  ) 

K3= ' II K > 

Oi  =  D I  K  I 
CGDS--CGD1  K  ) 

SGDS  *SGD< K  ) 

DO  200  L*l .NM 
NDl  *  HIM  L  ) 
l  A  =  I A  t  L  > 

L  6  I  A; L  ) 

Pl  * d ( i  ) 

'CD.  SGDI L  ) 

Nil  =■:,) 

DO  300  I  I  >  1  .  ND'< 

Fig.  A-l,  Subroutine  SGANTC 
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1  -MD  l  F  .  I  I  ! 

MM- (  1-1  )*N-(  1*1-1  >/2 
F  1  =  1  . 

1F.K3.E0. 1 2(  I  >  >U0  TO  36 
IF ( KB .EQ. I l(  I)  >F  I  =  ~  1  . 
15=1 

CO  TO  40 

36  IE( KA.EU. I3(  I  )  >  F 1 « -  l  . 

i  s  =  c 

40  DO  200  J0  =  1 . NDl 
■J  =MD  t  L  ,0J  ) 

MMM=MM+ J 

IF!  I  .  GT  .  J  )GO  TO  200 
F  J  =  1  . 


iF<  Lb.EQ. 12(0  )  ) GO  TO  46 
IF( L3 .EQ.  I  1(0  )  1FJ--1  . 

JS  = ! 

GO  '0  50 

46  IF(LA.£Q.I3<0> >FJ  =  -1  . 

FS  =  2 

50  IF  (  NIL  .  NE  . 0  >G0  TO  168 
N  I  L  -  1 

IFIK.fc'Q.l  )G0  TO  120 

I  N 0  = i l A - KA  ) * < LB-KA  )*( LA-K8  )*( LB-KB  ) 

IF' INP.EQ.01GO  TO  80 
C  SEGMENTS  K  AND  L  SHARE  NO  POINTS 

CALL  GGSCYL ( X<  KA  > , Y<  KA  > ,Z(KA  > ,X<  kB  > , V<  KB  > ,Z<  K8  )  ,X<  LA) , Y( LA  > ,Z<  LA  > , 
2  XU.  8 ) ,Y<  LB > ,Z< LB  > , RAOM , AM.DK ,CGDS , SGDS , DL . SGDT , I  NT. ETA, GAM 
3,n  1,1  ),P(  1 , 2  > .  P  (  2  ,  1  )  .  P  <  2 , 2  )  ) 

GO  TO  168 

C  SEGMENTS  K  AND  L  SHARE  ONE  POINT  (THEY  INTERSECT) 

80  KG  =  0 

J  M  =  F.  8 
JC  =  KA 
KF  =  ’. 

[NO  U  KB-  LA  >«<  KB-LB ) 

I F (  1ND . NE . 0  )G0  TO  82 
■)C  =  K3 
K  F  =  -  i 

■jm=f:a 
KG  =  3 

82  LG  =  3 
OP-lA 
L  F  -  l 

IF( L3.EQ. JC  >G0  TO  83 
OP  =  L  3 
t.  F  =  1 
LG=0 

83  SGN  =  KF  *LF 

CP  SI =( ( X(OP  )-X( JC  )  >*<  X< JM)-X( JC ) >♦< Y( JP )-Y( JC>  )•( Y(OM)-Y<  JC  )  ) 

2*<Z<  JP  >-Z<  JC  >  >MZ(0M)-Z<OC  >  >  >/( DK*DL  > 

CALL  GGMM(  .0,DK. .0,DL , AM , CGDS , SGDS , SGDT . CPS  I .ETA, GAM 
2.C(  1  ,  1  )  ,Q(  I  ,2 ) ,Q(  2, 1 > ,0( 2,2  >  > 

DO  93  KK=  1  .2 
KP  =  I AES < KK-KG  ) 

00  93  LL  =  1  ,2 
L  P  =  1 ABS ( L L -LG  > 

P  .  KF  .  LP  )  =  SGN *Q ‘ KK  ,  L  L  ) 

98  CONTINUE 
CO  10  168 

C  K-L  (SELF  REACTION  OF  SEGMENT  K> 

120  <H  1  '  .0.  .0) 


Pig.  A-l.  (Continued) 
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Cl)  2  =  (  .  0 . .  0  > 

1 F  <  CMM  .LE.tf.FGO  TO  150 
00= GAM* OK 
ZG=ZH/<  SG0S**2  1 
01 1 =ZG*< SGOS*CGOS-GO  )/2 . 

ZG  =  Zrl/(  SGDS**2  1 
Q 1 1 =ZG* ( SGDS'CG OS-GO  > / 2 . 

Q12=ZG*<GD*CGDS-SGDS  )/2. 

150  I 3CK- I SC( K  ) 

P 1  1  =  !  .0,  .0) 

P  1  2=  i  .0, .0) 

IF  '  I  JiCK  .  EQ  .  0  >G0  TO  155 
IF . BI-1.LE  AMFGO  TO  155 

gall  DSHELL (AM, BM,0K,CGDS.SGDS,EP2,EP .ETA. GAM, PI  1 . P 12  > 

t 55  ,1 1 =F i 1 +Q1 1 
Ql 2  -  P  1  2  ♦  Q 1  2 


Call  GGMM<  .0,OK, .0,DK,AM,CGDS,SGDS,SGDS,1 . 
2 . E  rA , GAM , P 1 1 .P12.P21.P22) 


01 

i 

1 

i  +Q  1  1 

-Jl 

*> 

_  -"l 

1 

2  +  Q  1 2 

p  • 

i 

.  i 

> 

*01  I 

Pt 

l 

,2 

=  Q  1  2 

?* 

> 

,  1 

) 

*012 

2 

> 

> 

=  01  1 

IF>  KA. ME .LA >G0  TO  160 
GO  TO  168 
160  P<  1  .  1  >  -  - Q 1 2 
PI  1 ,1  )*-Ql 1 
p  •  i  ,  M*-Q1  1 
P<  2.? >  =  -012 

168  C '  I  ,  J >  -C  <  I,0)*FI*FO*P(  IS.JSF 
IF  !  I  .ME  .0  >C<  0 , 1  )«C(  1  .0  > 

200  CONTINUE 

C  ADO  THE  LUMPED  IMPEDANCES  TO  THE  C  MATRIX. 

i see- 1 

301  lFU’LDl  ISEG).Mt.(0.  .0.  )>  GO  TO  302 
I F (  I 3EG . EG . I 2NM  >  GO  TO  350 

303  I  SEi-  - 1 6EG  *  1 
GO  10  301 

302  I  MODE  =  1 

307  IF'ISEG.GT.INM/  GO  TO  304 

IF  <  I2<  I MODE  ).EQ. I A  <  ISEGF  )  GO  TO  305 
306  I  MODE  -  I  MODE  *  1 

1 F (  I MODE . G  r . N  )  GO  TO  360 
GO  TO  307 

304  1eEG*ISEG  -  INM 

I  ■■  >  I2<  IMODE  >  .EQ  .  IB<  ISEG  >  )  GO  TO  305 
1 S E G -  I  o t G  *  INM 
GO  TO  306 

305  C  <  1 M DDF , 1  MODE  ) *C( IMODE , IMODE  )  ♦  ZLD(ISEG) 

GO  TO  303 

350  COwI 1  HUE 
RETURN 

360  .70  112(6.601)  1GEG ,ZLD( ISEG > 

601  FORMAT!  IX. ‘FATAL  ERROR:  UNABLE  TO  FIND  DIPOLE  MODE  TERMINALS*/ 
1  »  r OR  LUMPED  IMPEDANCE  ZL0(*,I3.*)  =  *,2E15.6) 

N  =  1 

..<)  TO  350 

r  *iD 


Fig.  A-l.  (Continued) 
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oooonon  nooo 


SUBROUT I N£  GGSCYL!  XA . VA.ZA. KB, YB.ZB.X1 , VI ,Z1 , X2 , V 2 , Z2 , RADM . AM 
2 . DS.CGDS ,SG0S . OT , SGDT , i NT , E TA . GAM , ° 11 .P12.P21.P22) 

COMPLEX  P11.P12.P21 .P22.CJA.EjB.EJ1 , E J 2 . ETA . GAM , C 1 .C2.CST 

COMPLEX  EGO  . CGDS.SCDS.SG0T.ER1 .ER2.ET1.ET2 

COMPLEX  ER11 .ER12.ER21 ,ER22 

COMPLEX  ET1 1 .ET12.ET21  ,£T22 

LOGICAL  AVIS, BV1S, VIS 

DATA  FP/12. 56637/ 

CA‘ i X2-X1 > /DT 
C  8  =  (  Y2-Y1  )  /DT 
CG  =  ( 22-Z1  )  /  DT 
CAS  =  - XB-XA  )/DS 
CBS-*!  YB-YA  )/DS 
CGS=( Z8-ZA  )/DS 
CC=CA*CAS+CB*CBS+CG*CGS 
I F ( A3S ! CC  )  .GT. .697  >GO  TO  200 
20  S2=( (1 -XA)*CASM Yl-''A>*CSS  +  iZl-ZA)*CGS 

THIS  VERSION  OF  GGS  WILL  TEST  IF  CYLINDER  SURFACE 
PREVENTS  "VISIBILITY"  BETWEEN  THE  TWO  SEGENTS. 


I F  (  I  NT . GT .0  )  GO  TO  301 
CALL  VISIBLE! XA.7A.Xt ,Y1 .RADM.V1S) 

IF!  .NOT. VIS )  GO  TO  301 
.'ALl  VISIBLE!  XB.  YB  .  X2  .  Y2  ,  RADM .  V I  S  > 

If  (  .NOT. VIS  )  GO  TO  301 
GO  TO  300 
301  INS-ZMINT/2) 

IF!  INS.LT.2  >INS  =  2 
I P  =  J  NS+ 1 
DE  LT-DT /INS 
T=  .  0 

D  SZ=*CC  *  DE  L  T 

PU0  .0.  .0) 

P 1 2s ( .0. .0) 

P  2 1  * !  .0. .0) 

3  22s  <  .0. .0) 

AMS*AM*AM 
S  G  N  I  . 

DO  100  I  N* 1 , I P 

OBSERVATION  POINT  (X  4  V  VALUES) 

T  X  ■=  Y  1  *  T*CA 

T  V  ~  Y 1  +  T  *CB 

X  &  V  VALUES  OF  THE  TWO  SOURCE  POINTS  ARE  XA  S.  YA  AND  XB  &  YB. 

DETERMINE  IF  THE  LINE  BETWEEN  OBSERVATION  POINT  AND  "LEFT"  END  OF 
THE  SOURCE  SEGMENT  INTERSECTS  THE  CYLINDER. 

CAlI.  VISIBLE!  XA  ,  YA.TX  ,TY  ,  RADM,  AVIS  ) 

C 

C  DETERMINE  IF  THE  LINE  BETWEEN  OBSERVATION  POINT  AND  “RIGHT"  END 
C  OF  THE  SOURCE  SEGMENT  INTERSECTS  THE  CYLINDER. 

CALL  VISIBLE! XB  .YB.TX.TY  ,  RADM , BVI S  > 

ET1 1 '0. 
t T21-0. 

£ T12-0. 
t  I  22=0. 

1 F ( AVIS .OR . BVIS  )  GO  TO  701 
GO  TO  750 

C  OBSERVATION  POINT  CAN  "SEE"  A3  LEAST  ONE  OF  THE  TWO  SOURCE  SEGMENT 
C  FNO  FOINIS. 


Fig.  A-2,  Subroutine  GGSCYL 


70  1  27  1  -  3  7 

7  lib 

»  >  :  >.  I*  !  *CA  -XA -SZ*CAS 
V'/  y 1 *T*CB-YA -SZ*CBS 

--z  i  .t*cg-za  -sz*cgs 
•»s.  x,Z**2WYZ** 2«ZZZ**2 

FAC  *  .  ■* 

If  ‘  Bi.iiT.  AMS  If  AC-C  CA»XVZfCB*YYZ  +  CG*ZZZ  l/RS 

:i  Ay  I".)  GO  TO  702 

GO  '  I  703 

ot'i:  P  v'A  r  I  !>N  'OINT  CAN  "SEE"  END  POINT  A. 

7  0?  -  '  VjR  ,  l  RS  *  Z7i  *»2  ) 

‘A-  I  ■  i  <  uAM*P 1  I 
•  i  .'API 

tW:»=L  )A-0^DS  ♦  ZZ 1  * E J  1  *CGDS 
t  T  :  I  =C(.*(  -E0  1*CGDS  )  ♦  E  AC*  E  R 1  1 
t  K  ? i  =  ZZl’EGl 
l  -  i  M  G*E J  1  »  <  AC*£R21 
>  8  «  ! S  )  GO  TO  703 

GO  I j  / 50 

UbD.  T  ;  ON  POINT  CAN  "SEE'  END  POINT  B. 

/»*  .>  T  <  R  S  ♦  ZZZ  *  *  2  I 

,  Sp(  - GAM*  ? 2 > 


f.  12=  1  G  B  /  R  2 
E  R  !  2  =  -ZZ2*EJ2 
El  1  ■■  -l.C  *  E  02  ♦  "AC.  *  E  R  1  2 
E R22 =■ E0B*SGP J  *  ZZ2-EJ2*CGDS 
ET22-CC*< -EJ2*;GDS'  *  FAC*ER22 
750  E  T l » E  T  1  1  ♦  E  T  1  2 
E  72  *  t r2 I  ♦  ET22 
C=  i. »  SGN 

IfCIN.EQ.l  .OR.  IN. EG.  I?  1C=1  . 

E  GP =  C£  XP ( GAM* ( 9 T - T >  > 

1  --C.  *(  EGD-l  .  /EGO  )/2  . 
t  GO  =  CE  VP ( GAM*T  ) 

'",=C*t  EGD-l  ./EGO)'?. 

P 1 1 1 1+FT1*C 1 
P  I  2  =  MM  T  1  *C2 
P21-"32l-"ET2*Cl 
P22=  322*ET2*f.2 
r=T  f  )E  LT 
vz=sz*csz 

100  G  ON  -  -GGN 

CST=-ETA"DELT/t  3. *FP*SGDS»SGDT  > 

P 1 1«CST*P  1  1 
P|2=CST*P12 
F ?1=CST*P21 
P22=CST*P22 
RE  To  IN 

200  SZ1  =  .  v  1  -VA  > *CAS  +  <  VI  -VA  > *C8 S *•  <  ZI  -ZA  )*CGS 

PHI  =  3'jRTC ( X 1 -XA-SZ1*CAS  ) **2 * ( Y 1  - YA-SZ 1 *C B S  > * * 2 ♦ <  Z 1 -ZA-SZ 1 *CG S  )* 

3Z?  =  SZ! ♦ D  T  *  C  C 

R  H2 “ SORT ( ( X2-XA-SZ2*CAS  )  *  * Z * ( Y 2 - YA- S Z2 *C BS  > * * 2 ♦ ( Z2-ZA-SZ2 *CC S  )* 
D Li P  =  :  RH1+RH2  I/7. 

in  ODD .GT . 20. *AM  .AND.  INT.GT.01GO  TO  20 
I F  : DID  .  L  T . AM  I D 3 0  *AM 

CALI  CGMMC  .  0 , 0 3 . S Z 1  .  SZ2 , ODD  ,CGDS  .  SGDS , SGDT , 1  . 

2  .  I  TA  .GAM, P 1 1  .  P  1  2  .  P21  ,T 22  1 
KF  T L «  V 

300  S'  --GORTC  1  .  -CC*CC  ) 

<"  AD  *  (  i.GS*l.B  -  LBS  *rc.T  )  /GO 

Fig.  A-2.  (Continued) 
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CBD»(OS*rG-CG;*CA)/SS 
-GD  =  t  CBS'CA -CAS*CB ; /Sf 

0K°A9S(DV  )*CAPMYi"YA,,'C8D  +  <Z1':2A>*CGD 


If ( OK . LT . AM  )  DP  =AM 
w'=XA*SZ*CAS 
'■  2  =  /A»S2*CBS 

22=za+sz*cgs 


VP  1  *■  K  1  -OK'CAO 
VP  1  --V  1  -  OK *C BO 
’I  I  -  21  -  DK*CGD 
>«  ’•-GBS»CGD-CGS*CBO 
'-•BP»CGS«CAD-GAi*CGD 
«  GP»CAS«CBO-CBS*CAD 
P I “GAP*( XP  1 -*Z)+CBP*( 
T I  P  1  /  SS 
'II  ■  1  I  *CC-SZ 

GGMMISI  ,SI*DS  U 
2  PI1.P12.P21  P  '  2  > 
PtTURN 
:  no 


VP  1 -YZ >+CGP*( zp 1 -zz  ) 

.Tl+DT.DK, CGDS.SGDS.SGDT.CC, ETA , GAM 


Fig.  A-2,  (Continued.) 
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include  the  x  and  y  end-points  of  the  source  and  observation  segments, 
as  well  a:;  the  radius  of  the  cylinder  in  meters.  The  z  values  of  the 
end-points  ore  not  required  since  the  axis  of  the  cylinder  is  assumed  to 
lie  on  the  z  coordinate  axis,  as  for  the  cylinder  reflection  subroutines 
written  by  Skelman.  VISIBLE  will  return  the  logic  variable  VIS  which 
will  be  "FALSE"  if  the  surface  of  the  cylinder  intersects  a  straight 
line  drawn  between  the  source  point  and  the  observation  point.  As  with 
CSS,  line  16  will  test  if  the  source  and  observation  segments  are 
parallel  to  each  other  with  small  displacement.  If  so  GGSCYL  calls 
GGMM  just  as  GG3  does.  If  closed-form  impedance  calculations  have  been 
selected  by  the  User,  which  is  done  by  setting  variable  INT  equal  to 
zero,  statements  23  thru  26  will  determine  if  the  cylinder  surface 
obstructs  direct  radiation  between  the  two  segments.  If  so,  the  use  of 
subroutine  GGMM  is  not  permitted.  If  the  mutual  impedance  is  to  be 
calculated  via  Simpson's  rule,  GGSCYL  will  again  determine  if  the 
surface  of  the  cylinder  will  interfere.  At  this  point,  GGSCYL  imple¬ 
ments  the  concept  introduced  by  Ekelman  in  the  reflection  subroutines 
SG  where  the  source  segment  current  is  considered  to  be  concentrated  at 
points  on  either  end  of  the  source  segment  as  shown  in  Fig.  A-3.  The 
end  points  of  the  source  segment  are  labeled  A  and  3  for  convenience. 
Subroutine  GGSCYL  then  determines  the  field  which  arrives  at  various 
points  T  that  lie  on  the  observation  segment  from  the  two  end  uoints  of 
the  source  segment,  provided  that  the  cylinder  surface  does  not  intersect 
the  linn  between  the  source  segment  end  point,  A  or  B,  and  the  observa¬ 
tion  Dolnt,  T.  If  it  does,  the  field  contribution  from  that  end  point 
is  zero.  GGSCYL  then  integrates  the  r;. eived  field  along  the  observation 
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a  <  XB, YB, ZB  ) 


-vi€.  A-3. 


interaction  between  source  and  observation 
segments  in  subroutine  GGSCYL 
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segment  just  as  the  original  GGS  subroutine  does  to  determine  the 
mutual  impedance  term  for  that  pair  of  segments. 

Subroutine  VISI3LE,  listed  in  Fig,  A -4,  is  frequently  called  by 
other  subroutines  to  determine  if  the  surface  of  a  circular  cylinder  of 
radius  RAD,  whose  axis  is  coincident  with  the  z  axis,  intersects  a 
straight  line  drawn  between  points  P  and  R.  If  the  line  does  not 
intersect  the  surface  of  the  cylinder,  the  logical  variable,  VIS,  is  set 
to  a  value  of  "TRUE".  Since  the  problem  is  independent  of  the  value  of 
z,  the  cylinder  is  reduced  to  a  circle  in  the  z  =  zq  plane.  Because  of 
its  repeated  use  by  many  other  subroutine  programs,  VISIBLE  was  designed 
to  execute  as  fast  as  possible.  The  geometry  of  the  problem  is 
illustrated  in  Fig.  A-5.  First  the  distance  of  points  P  and  R  from  the 
axis  of  the  cylinder,  j OP  |  and  | OR  f ,  is  determined  and  if  they  are  less 
than  the  radius,  at  least  one  of  the  points  lie  inside  the  cylinder. 

This  is  possible  for  the  monopole  antenna  model  used.  In  that  case,  the 
points  can  not  "see"  each  other  and  the  value  of  VIS  Is  set  to  "FALSE", 
Next,  the  direction  cosines  of  the  line  thru  P  andRare  determined. 

Point  Q ,  which  lies  on  line  PR,  Is  determined  such  that  the  length  of 
line  |  OCJ  |  represents  the  minimum  distance  between  the  cylinder  axis  and 
the  line  PR.  This  requires  that  line  OQ  be  perpendicular  to  line  PR. 
Using  tliis  fact  and  a  theorem  from  analytical  geometry  [12]  which 
determines  the  distance  between  a  plane  (or  line,  sketched  as  a  dashed 
line  in  Fig.  A-5) ,  normal  to  a  line  (PR)  that  intersects  the  plane  at  a 
point  (?) ,  and  the  origin,  the  distance  PQ  can  be  very  quickly  deter¬ 
mined.  This  allows  one  to  locate  point  Q.  Using  the  relationship  of  a 
right  triangle  (AoPQ)  the  minimum  distance  between  the  line  FR  and  the 
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S 1 1 K  K '  • :  i  I  INF  VIS:  F  1  t  <  P  X  .  I  V  ,  R  X  .  R  V  ,  HAD  .  V I  S  > 

C 

C  I  H  I  '  Ml  OUT  1  N  E  DETERMINES  IF  A  STRAIGHT  LINE  DRAWN  THRU  JOINTS 
C  R  AND  I;  WILL  INTERSFU  THL  CIRCLE  WITH  RADIUS  A  WHOSE  CENTER  IS 
C  AT  THE  ORIGIN.  IF  THE  LINE  PR  DOES  NOT  INTERSECT  THE  CIRCLE, 


C  THEN  VI:'."  «  “TRO:". 

C 

I  '  C  I C  A  L  v  '  S 
?  I  AL  P  ?  ,  I  R  2 
r  ^ -  V,  ♦  P  1  ’  P  7 
lR2-t>*rx  *  r, ,  *  r  v 
R AlV  -RAD*  RAD  -  .  I 

IEiLM:  LEVADA  .OR  LR2.lt. RAD2  )  GO  TO  900 
US=  .1  RUE  . 

C  COMPUTE  DIRECTION  COSINES  THRU  P  4  R. 

D  L ' G  0  F  T  <  <  R  > - P  V  ■  *  *  J  ♦  (  R  V  - P  V  )**2  ) 

I  E  (  D.  r  n . 0 .  >  RE  TURN 

•' :  *  t  >  .<  f  x  -  -  la 

.  2  *  l  <  1  -  r  Y  )  '  in 

C  COMPILE  DIRECTION  LENCIH  Or  PERPENDICULAR  LEG  Sf>  ON  LINE  °R. 

C  (SIMI.AR  TO  SZ  CALCULATION  IN  SBR  GGS .  ) 

SP  *  -  :  <»0  1  -  P V  *C2 

C  FIND  LOCATION  OF  °EllPENBICU_AR  INTERSECTION  Q 
I  X =  U  *  SP  *  P  X 

C  TEST  IF  0  lies  BETWEEN  P  AND  R. 

IF  (  ’  < . GT  OX  .AND.  QX.GT.RX  .OR. 

1  RX.GT.QX  .AND.  OX.GT.PX)  GO  TO  100 

o  /  •=  c :  *  s  p  '  p  v 

>  PV  GT.OV  .AID.  OV.GT.RY  .OR. 

1  P/.GT.GV  .AND.  QY.GT.PV)  GO  TO  100 
C  CLOSEST  POINT  TO  CYLINDER  DOES  NOT  LIE  BETWEEN  9  AND  Q. 

PS  TURN 

100  S'’2'bT  *SP 

C  COM° U T E  LENGTH  OR  F E RP E ND i C JLAR  LINE  FROM  LINE  D»  TO  CYLINDER  AXIS. 
DP  *  SOP  T ( LP2-SP 1  ) 

C  DETERMINE  IE  INTERSECTION  WITH  CYLINDER  OCCUP.ES. 

1  F ( D J  .  l  E  . PAD  )  V S * . FALSE . 

RETURN 

C  PUiNT  IS  INSiDE  THE  CYIINDER;  OK  IE  D ART  OP  ’MONOPOLE"  ANTENNA 
C  HOWEVER,  POINT  MUST  NOT  BE  “VISIBLE"  TO  OUTSIDE  POINTS. 

900  VII-  .FALSE  . 

P  E  T  U  ?  N 
END 


Fig.  A~4.  Subroutine  VISIBLE 
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R* 


Fig.  A -5.  Geometry  concepts  of  subroutine  VISI3LS 
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cylinder  axis  is  rapidly  computed  and  compared  to  the  radius  of  the 
cylinder.  Thus,  if  the  distance  |  OQ  }  is  less  than  the  radius  of  the 
cylinder,  line  PR  intersects  the  cylinder.  There  is  a  situation,  however, 
where  both  points  may  be  located  on  the  same  side  of  the  cylinder,  as 
for  points  P*  and  R'  shown  in  Fig.  A-5.  In  this  case,  although  line 
P'R’  intersects  the  cylinder,  the  cylinder  surface  does  not  intersect 
the  line  segment  between  points  ?' and  R'  and  therefore,  ?'  and  R'  can 
still  "see1'  each  other.  However  for  this  case,  point  Q',  which  marks 
where  line  P'Q'  is  closest  to  the  origin,  does  not  lie  between  the 
points  P'  and  R’.  Thus,  in  the  subroutine,  if  point  Q'  does  not  lie 
between  the  two  points,  and  both  points  are  located  outside  the  cylinder 
(tested  at  the  beginning  of  VISIBLE),  then  the  two  points  will  be  visible 
to  each  other,  regardless  of  how  close  the  line  thru  the  two  points  is 
to  the  cylinder  axis.  This  is  the  purpose  of  the  test  at  statement  29 
of  the  subroutine. 

Subroutine  SG  determines  the  delta  impedance  matrix  terms  due  to 
reflections  from  the  surface  of  the  cylinder.  This  subroutine,  listed 
in  Fig.  A -6,  is  the  routine  written  by  Ekelman  except  that  a  call  to 
subroutine  VISIBLE  Is  made  at  statement  159.  If  the  surface  of  the 
cylinder  is  determined  to  intersect  a  straight  line  between  the  source 
segment  end-point  and  the  observation  point,  then  there  is  no  field 
which  arrives  at  the  observation  point  due  to  reflections  from  the 
cylinder  surface. 

Subroutine  SGSURF,  listed  in  Fig.  A -7,  is  the  subroutine  which 
calculates  the  delta  Impedance  matrix  terms  due  to  boundary-layer 
surface  ’wave  diffraction  from  the  surface  of  the  cylinder.  The  first 
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.US'  ji)  r  1  nf  sc,:: i  ico  ,  ia,  ib .  inm,  int .  isc .  1 1 . 12. 1 3.  ja.  jb  .md . n .nd  . nm . 

2!-  ‘  .Am.BM.^ .CGO.LMM.D.EP2.EP3.ETA  ,  F  HZ  ,  GAM  .  SGD  .  X  .  V  .  7  .  Zl  D.2S.A.B.WL  ) 
C  *  *  St,*.  IS  THE  SAME  AS  SG  EXCEPT  THAT  SUBROUTINE  “VISIBLE"  IS  CA-LtD 

*  -  '  -  r-ECLUJE  cOjPlING  ANALYSIS  BETWEEN  TWO  SEGMENT  ’UNI'S  HIDDEN 

*  *  FRiM  EACH  OTHER  BY  THE  CYLINDER. 

COMP-  EX  ZG ,  Zw . ZS .EGD.GD.CGDS .SGDS.SGDT.B01 

COMP L E X  Pll.P12.P21.P22.011.Q12.Q21 , Q22 . E P2 . E P . ETA . GAM , E P 3 

complex  epsila.cwea.8Eta.zarg 

COMPLEX  P(2,2).Q(2.2)  ,CGD(  1  ) ,  SGD (  1  ) .  C (  I C J . I C 0 > . ZL D<  1  ) 

COMPLEX  EX1 ,EY1 ,EZ1 ,£X2,£Y2,EZ2,ET?R,ETPP,PH,£RPR.ERPP 
COMP-EX  EX.EY,EZ,ERX.ERV,ERZ,ET1,ET2.EIPR.EI?P.C1.C2.CST 
LOGICAL  VIS 

DIMENSION  UN< 2  )  . UB( 2  )  . VI : 3 > 

DIMENSION  X(  I  )  .  Y<  1  ) ,Z<  1  > ,Df  1  )  . I  A;  1  I  ,  I  B (  I  )  . MD 1  I NM . 4  ) 

Dimension  i i ( i  >,1211 > . 1 3 < 1  > . oai 1  > , j b < 1 > . n d < 1  > ,  1  sc<  1  1 
DATA  E0.TP ,U0/8 -85  4E-I2 .6.28318 . 1  .2566E-6/ 

DATA  FP/  12. 56637/ 

2  FORMAT! 3X . *AM  =  * . E 10 . 3 . 3X . *DMAX  =  * . E 1 0 . 3 . 3 X . * DM 1 N  =  *.E10.3> 


*  WARNING: 

*  SG  IS  THE  ONLY  SUBROUTINE  WHERE  THE  VA.UE  OF 

*  THF  CYLINDER  RADIUS  "A"  IS  IN  UNITS  OF  WAVE 

*  LENGTHS.  NO'  METERS 

* 


E  P  *  E  R  3 

!CC»'N*N  +  N  )  / 2 
DO  I  1  0  =  1  .N 
DO  1 J  I =1  ,N 

10  C (  I . 0  )  =  ( 0 .  , 0 .  ) 

11  CONTINUE 
ZS“<  .0, .0) 

I  F 1 CMM . L  E .0.  >G0  TO  12 
OMEGA^TP *F  HZ 

E  P  S  I  LA  =  CMPL  X ( E0 . -CMM* 1 . E  6 /OMEGA » 

CWEA=<  .0,  1  .  >«CMEGA*EPSKA 

bETA  =  OMEGA*SQFT<  U0)*CSORT< EPSILA-EP  ) 

ZARG=BETA*AM 

CALL  C B E S (  ZARl.  ,801  ) 

7  G - B  F  TA*  B0 1 / CWEA 
12  ZH " ZS / <  TP  * AM* GAM  ) 

DMI N= I . E  30 
DMAX -  . 0 
DO  20  J=1  . NM 
K  I A  (  0  ) 

L  =  I  B  (  J  ) 

0(0)  * SQRT  <  (X!K)-X!L)>**Z  +  (V(K>-Y(L>)**2  +  (Z!K)-Z(L!>  =  =  2> 
1 F ( Oi 0  )  .LT.DMIN  >UMIN=D(0  > 

IF ! D( 0  > . GT . DMAX  )DMAX  =  D( 0 > 

EGD=CEXP! GAM* D ( 0  >  > 

CG0!0  >  =  < EGD+  l  .  /EGO  1/2. 

20  SCO! 0 >  =  ( EGO- I  . /EGD >/2  . 

IFiDMIN.LT.2. *AM  )GO  TO  25 
I F ( CABS! GAM*AM ) . GT . 0.06 )GO  TO  25 
I F ( CAES! GAM*DMAX  )  GT.3.  )GO  TO  25 
I F ! AM . GT.0. )GO  TO  30 
25  N  ~0 

UR  1  T.  1  6  .  •  'AM  .  DMAX  .  P.MIN 
RETURN 

30  )u  200  K-l.NM 
NO'  —  'If* !  K  ' 


Fig.  A -6.  Subroutine  SG 
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VA -  ' A<  K  I 
HE= I3>  K  > 

DK  =  D( K  > 

CGi)3  *CGD<  K  ) 

sens  «SGD< K  ) 

00  j  fi-  L  =  1  ,  NM 
NDL  -  ND(  L  ) 

L A= I A i L ) 

I.  B  =  I  3  (  L  ) 

0  L  =  D  (  L  ) 

SGOT  -  GGD(  L  ) 

N  I  L  =0 

DO  2 00  11=1 , NDK 
I  ~MI  (  «  ,  l  I  ) 

MM=  1-1  >*N-(  1*1-1  )/2 
F  I  --  1  . 

IF(  F a .EO.  I2<  ]  > >G0  TO  36 
1  F  (  1, 3  .  E  Q .  1 1  (  1  >  )  F  I  =  -  1  . 

I  S=  1 

GO  TO  40 

36  1 F ( KA . EQ .  13(  I  )  > F  1  =  —  1  . 

13  =  0 

40  DO  2  30  0  J  =  1  ,  N  D  L 
J=H')(  l  ,0  J  ) 

MMM=HM+0 

C  IF (  1 .GT . }  )G0  T0  200 

F  0  =  1  . 

IrL3.EQ.i2<J;  >G0  TO  46 
IF  <  LB.tQ.  11(0  >  >FO  =  -I  . 
JS=i 

GO  TO  50 

46  1 F ( LA . EQ .  I 3<  J  )  >F J  =  -l  . 

J  S  *  2 

60  IF ( NIL . NE ,0!GO  TO  168 
N  I  L  =  1 
CC 

CC  BEGIN  EQUIVALENT  TO  GGS 
CC 

CA=  :  ><(  LB  >~X<  LA  )  )/DL 
CS“ :  y i LB  ;  -Y( L A  >  )/DL 
G - <  Zi  LB  :-Z<  LA  )  )/DL 
I  NS  =  ? * (  I  NT/ 2  ) 

IFF  I nS . LT . 2  > INS=2 
1  P  •  1  N  $  *  1 
L't  I  T  =  DL  /  I  NS 
T 

p  n  =  (0.  ,0.  i 
f  l :; = ;  0 . ,  0 .  > 

F  2  i  =  i  0  .  ,  0  .  ) 

P 22=10. ,0. ) 

P>  1  •  1  )  =  <0. ,0,  ) 

Ft  1  .2 >*(0.  .0.  > 
p<  2. 1 )  =  '0.  ,0.  ) 

F  1  2 . 2  >  =  (  0  .  ,  0 .  > 

SGN= - 1 . 

DO  100  I N= 1  ,  IP 
XC  =  (  XI LA ) -T*CA )/WL 
VC  =  (  y(l  A  >+T-CB >/Wl 
ZC  =  ■  2(  lA  )  +  T*CG  )  /WL 
E  T  :  =  I  0  .  ,  0  .  ) 
t  r .'  •=  ( 0 . ,  0 .  > 

DO  103  JEP-  1  .  2 
I ;  I  P  .  F.O  .  2  >G0  10  8  3 
v  SE  =  x  (  KA  5  /Wl 


Fig.  A-6.  (Continued) 
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V'jtt  W  t.A  >/WL 
kA  '/WL 
GO  TO  a 4 


83  -  U-  >(  ,:ij  •  WL 
*  .P  •  y  :  f  B  '  /Wi 
0  0  ■•fll/VT 

84  i  Of,  i  i  NIJE 

c 

C  VISIBILITY  TE  S i  BETWEEN  SOURCE  SEGMENT  END  POINT  AND  OBSERVATION  PT. 
C 

CALL  VISIBLE  <  XSR.YSR.XC.YC.A.VIS  > 

II  <  .  kOf . VIS  )  GO  TO  102 
CC 

CC  RE  ’  TC  FINDS  THE  REEL  PT  ON  THE  CYLINDER 
CC 

C  A  _  L  Ri N  P  T  C ( XSR . YSR.ZSR, XC. YC.ZC.A.B , XR, YR ,ZR , VR. UR. VI . VIM  > 

CC 

CC  GN*  FINDS  THE  F I E  L  D  AT  THE  RE  F  L  PT  WITH  K-SEG  AS  SOURCE 

CC 

XRR-  *.R  *UL 
Y  P  R -  FR*WL 
2RR=/R*WL 

CAl L  GNNF ( X( KA  ) , Y<  KA),Z( KA>.X( KB  ) , Y<  KB  >.Z( KB  > ,XRR, YRR . 

2  ZRR.AM.DK.CGDS.SGOS.ETA.GAM.EXl ,EY1 ,EZ1 .EX2.EY2.EZ 2. JEP  > 

CC 

CC  NANDB  FINDS  THE  NORMAL , B 1  NORMAL .AND  TANGENT  VECTORS  AT  REEL  PT 

CC 

CALL  NANDBIA.B.UN.UB.VR) 

CC 

CC  BEGIN  TO  FIND  FIELD  AT  XC.YC.ZC  DUE  TO  REFLECTED  FIElD  AT  XR.YR.ZR 

CC 

P  '  =  7? >2. 

TF  I  =  TP 

C  FH  I  -UNI  1 '  *V1( 1 )»UNI 2  >*71  I  2  > 

SXN»XSR-XR 
SVN*  fSR-VR 
SZN  =  ESR -ZR 

SMAG-SQRTI SXN«SXN  +  SVN*SYN  +  SZN*SZN  > 

SXN» SXH/SMAG 
SV  N  =  SYN/SMAG 
SZN=SZN/ SMAG 

C  TESTING  FOR  B./0.  IN  AT AN 2 
TCMP=SXN*UB<  1  >*-SYN*UB(2) 

IF  l TEMP .EQ.0. 1  GO  TO  5 
WR  «ATAN2( TEMP , SZN ! 

GO  TO  6 

5  W'II->tF.0 

6  SW=  3  I N ( WR 1 
CW«COS(WR > 

SST2 - 3W*SW+CW*CV*CTH I *CTH I 

PH02-SMAG 

CSV»COS( VR ) 

SHV=3IN( VR  > 

OD»SORT( 8*B*CSV*CSV+A*A*SNV*SNV  > 

RG-DO*DD*DO/A/B 

R  HO  I =SMAG*RG*CTHI/( RG*CTHI *2. *SMAG*SST2 > 
i,IPRX  =  SIN<  WR  -P  I  /  2  .  t*UB<  I  ) 

U  >  PRY-SIN, WR-PI/2.  >*UB<  2  > 

U  I  PR  Z  -  COS  ( WR - P  I  / 2  .  ) 

UIPPX’SYN*U!PP2- SZN*U /PRY 

Pig.  A -6.  (Continued) 
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u  m'  r-s/N«u i  pp  <  -sx..  -n  i  pr_i 
U  IPP.''*IXN*U)N<  ■  '.""i  'M 
UP.l’P  <-U  IPRs  *V!  '  1  )  -  1 1  1  ■/  1  '  . 

ORPP»  =  iiiPR.i*vi:  i  )  1 1 -'R ) 

URPP  >UIPRX*V1  (  2  )-  llr>R  Y*M  ;  ) 

EX--E  <1 
EY  =  E  yi 


tr=Ezi 

DO  1.01  I  MODE  '1.2 

E1:'R-(U1PRX*EX<-UIPRY*CY  !;lRZ*EZ) 

E  )PP-(ijippx*EX*l'IPr'V*,r  ‘  ! '  "Z^EZ  ) 

PH*CiXP<  CMPLXt  .0.  .  -TP  I  *7  ’■  ■■  /  SORT!  <  RHOl  +  VIM  >*(  RH02+VIM  >  ) 

E KP R - -SORT ( RHOl *RHO:  ,  "'-i*C  1  p .0 
E  RPP “SORT ( RHOl *RH02  )  V  .  •  1  .  , 

ERX«ERPR*UIPRX«-ERPP*** 
rRV-ERPR‘UIPRY  +  ERPf-  *  0  -  ■< 

EPZ-ERPR*UIPRZ*-ERPP*UR!V2 

CC 

CC  DOT  THE  FIELD  WITH  THE  _-iIG  DIRECTION 
CC 

in  I  MODE . E  Q . 1  >ET1=ERX*CA+£RY*CB+ERZ*CG+ET1 
I F ( MODE .  E  Q  .  2  1 E  T  2  =  £  R  X  1  C  A  ♦  E  R  Y  *  C  B  +  E  R  Z  *  C  G  +  E  T  2 
EX  =  E  <2 
E  Y  =  E  Y  2 

101  EZ*EZ2 

102  CONTINUE 
CC 

CC  ASSIGN  INTEGRATION  WEIGHTS  AND  MODAL  CURRENT  VALUES 

CC 

cwi  =  :< .  *  SGN 

rniN.EQ.l  .OR.  I N . EQ . IP  >CWT » 1 . 

E  GD“CE  XP ( GAM* ! OL-T ) 1 
C  1  «cv:  *(  EGO-1  .  .'EGO  1/2  . 

EGD=CEXP<  GAM'T  1 
C2=CWT ‘ ( EGO-1 . /EGO ) ■ 2 . 

PU^l  l+ETl»Cl 
P  1Z  =  P  12*ET1*C2 
PZl-RZnETZ^Cl 
P  _2  -  5«1?+ET2*C2 
T-T-Of.LT 
100  SC.N--SGN 

CRT  -  -DELT/I 3 . *SGDT > 

P  1  1  =  CST*P l 1 
•»l?  =  CST»Pi2 
P  2 1 =CST*P21 
P22=CST*P22 
P  I  1  P  1  1 
P<  1  .  ?  >»P12 
P  (  2  ,  -  P  2  1 

P!  2.  M  =  P22 

168  C<  I  ,0  )  -  C  <  I  ,0  >*';!*FJ*P(  IS, OS) 

C  ir(i.NE.J)C(J,l)*C'[,J) 

200  CONTINUE 
RETURN 
END 


Fig,  A -6.  (Continued) 
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SUBROUTINE  SGSJRf i ICO , IA , I B , I NM . I  NT , I  SC , II , I  2 . i 3 , OA , OB , MD , N . ND , 

1  NM,AM,BM,C,CGD,CMM,D,EP3,ETA,FHZ, GAM ,SGD,X.V,Z,ZLD, 

2  ZS , RADM ,WL . NC , F  AC ) 

COMPLEX  CSGDS.EN1 ,EN2  ,  EMAGS , EMAGL 
COMPLEX  ZH,ZS,EGD,CGDS,SGDS,SGDT,B01 
COMPLEX  PI  1  ,P1 2.P21 ,P22 ,EP ,ETA,GAM,EP3 
COMPLEX  EPS1LA,CWEA,BETA,ZARG 

COMPLEX  P<  2 ,7  )  ,CGD(  1  )  ,SG0(  1  )  ,C(  ICJ, ICO  > ,ZLD(  1 > 

COMPLEX  ET1 , ET2 ,C 1 ,C2,CST 
LOGICAL  VIS 
DIMENSION  NCI  1  ) 

DIMENSION  X< 1 ) ,Y<  1),Z( 1 ),D( 1 > . IA< 1 >. IB( 1 >.MD! INM.4  ) 

DIMENSION  Il<  1  >,  12<  1  >,I3( 1>.0A( 1 >,OB< 1 >,ND<  1  >.ISC<  1  > 

DATA  E0.TP , U0/3. 8S4E- 12,6. 28318, 1 .2566E-6/ 

2  FORMAT! 3X , *AM  =  *  ,  E  10 . 3 . 3X , -DMAX  -  * , E 10 . 3 , 3X , *DMI N  »  *,E10.3> 

E  P  =  E  P  3 
DC  11  0  =  1, N 
DO  10  I  *  1 , N 

10  C(  I  .O >»<0. ,0.  > 

11  CONTINUE 
ZS*i .0. .0) 

If (CMM.LE.0. )GO  TO  32 
OMESA=TP«FHZ 

E  pS I LA-CMPLX( E0, -CMM* 1 . E 6/OMEGA  1 
C WE  A  =  (  .0, 1 . >*OMEGA"EPSILA 


BETA -OMEGA *SQR I ( U0  >*C SORT( EPSILA-EP  > 

ZARG=BETA*AM 

CALL  CBES ( ZARG , B01 > 

ZS-BETA*B01/CWEA 

12  ZH=ZS/(TP*AM«GAM) 

DM I N  3 1 . E30 
DMAX=  .  0 
DO  20  0=1 ,NM 
K=IA( 0  ) 

L  3  I  B  !  0  ) 

D  (  J ) =SQRT  ((X(K)-X(L) )**2  +  ( Y( K ) — Y  < L  )  )**2+(Z(K)-Z(L  >  >*  *2  ) 
IF<D<0>.LT.DMIN>DMIN«D!U> 

I F ( D ( 0  )  .GT.DMAX  1DMAX-D10  ) 

EGD- C EXP! GAM*D( 0  )  ) 

CGD ( J  )-<  EGD  +  1 . /EGD  )/2  . 

20  SGD( 0  >“( EGD-1 . /EGD  )/2. 

IFIDMIN.IT. 2. *AM)GO  TO  25 
I F ( CABS( GAM* AM  >.GT.0.06)GO  TO  25 
I F ( CABS ( GAM*DMAX  ) . GT . 3 .  >GO  TO  25 
IFCAM.GT.0. )GO  TO  30 

25  N-0 

WRITE! 6,2  >AM , DMAX , DMI N 
RETURN 

30  DO  200  K3 1 , NM 
NDK-ND ( K  ) 

KA* I A( K  ) 

KB»  I  B ( K  ) 

DK>D< K  ) 

CGDS-CGD! K  ) 

SGDS-SGD! K  ) 

C3GPS--CGDS/SGDS 
DO  200  L  3 1 , NM 
NDL-NDl L  ) 

LA-IA1L  ) 

LB-  I3<  L  ) 

Di  -DILI 


Pig.  A -7.  Subroutine  SGSURF 


SGOT-SGOI L  ) 

NIL  =0 

DO  2.30  11  =  1  ,  NDK 
I =MD<  K  .  1 1  ) 

MM=<  1-1  )*N- ( 1*1-1  1/2 
F  I  =  1  . 

IF! K3.EQ. I2( I >  >GO  TO  36 
I F  < K3.EQ. 1 1(  I  >  >FI«-1  . 

ISM 

GO  TO  40 

36  IF; KA.EQ. 1 3  <  I  >  )F  I  =  -l  . 

IS*2 

40  DO  200  OJ  =  i  ,  NOL 
J= MD! L ,00  > 

I ?  I  *  l  2<  I  ) 

120  =12(0  ) 

IF ( NCI  i2l  )  .EQ.0  .OR.  NO ( 120)  .EQ.0)  GO  TO  200 
C  I F (  I . GT . 0  >G0  TO  200 

F  0  =  1  . 

I F  < L8.EQ. 12(0  )  )G0  TO  46 
IF ( L3.EQ. I K  0  > >F0  =  -1  . 

0C=  1 

GO  10  50 

46  I F  ( l  A .  E  Q . 13(0  >  )  F  0  =  - 1 . 

0  3  =  2 

50  I F ( NIL . NE .0  )GO  TO  168 
N  I  L  M 
CC 


CC  BEGIN  EQUIVALENT  TO  GGS 
CC 

C  COMPUTE  DIRECTION  COSINES  OF  EXPANSION  MODE  SEGMENT  L 
CA  =  <  X ( LB  )-X( LA  > >/DL 
CBM V( LB >  —  V < LA)  >/0L 
COM  Z( LB  )-Z( LA)  >/DL 

C  COMMUTE  DIRECTION  COSINES  OF  TEST  SEGMENT  X 
CAS  =  <  X<  KB  )-X( KA)  >/0K 
CBS*< V( KB  ) - Y  < KA  )  >/0K 
CC5  =  ( Z( KB  )  -  Z( FA  )  )/DK 
I  NT  1*1 NT 

! F (  NC( LA ) . EQ .0  .AND.  NC(KA).EQ.0  )  1NT1=20*INT 
I F (  NC< KB J.EQ.0  .AND.  NC(LB).EQ.0  >  INT1=20*INT 
I F (  NC< LA  > . EQ . 0  . A N 0 . N C ( <B ) .EQ.0  >  INT1=20*INT 
I F (  NC l KA  > . E  Q . 0  .AND.  NC(LB).EQ.0  )  INT1»20*INT 
I NS  =  2* ( INTI /2  ) 

I  F (  I  NS . IT . 2  > I N$“2 
I  P*  I NS+ 1 
DEL  T  =  DL / INS 
T  =0  . 

.0. ) 

Pi2=:0. ,0.  ) 

P21*«  0. .0.  > 

P  2  2  « <  0 .  ,  0 .  > 

P  (  1  .  i  )  ( 0 .  ,  0 .  > 

PM  .  2  >  =  <  0.  .0.  ) 

P  t  2  ,  l  >  *  (  0  .  ,  0  .  ) 

P ( 2 . 2  ) *( 0 . .0.  ) 

SGN=  -  1  . 

DO  100  I N  =  1  , IP 
PX  =  i  <!  l.A>‘T*CA) 

PY-t  Y( LA  )  *T*CB ) 

PZ-  ' 21  LA ) *7*CG ! 

C  LOCATE  SOURCE  CURRENT  AT  CENTER  OF  SOURCE  SEGMENT 
P P  X *  X ( KA  >  *  Dk  * . 5 *CA5 
P T  Y  ■=  f ’  (  k A  )  *  OK*  .  5  *C6 5 
PlZ-T<  KA  )  ♦  DK*  .  5 * C C S 
C 
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DO  NOT  COMPUTE  SURFACE  E-FIELD  IF 


C  DETERMINE  IF  P  CAN  “SEE*  PP . 

C  THE V  CAN  “SEE"  EACH  OTHER. 

C 

CAL L  VISIBLE*  PP X , P P Y . P X . P Y , RADM , V I S  > 

1 F  <  .NOT. VIS  )  GO  TO  261 
£ N 1 *CMPL  X <  0 . .0.  ) 

£N2*CMPLXI 0 . ,0.  ) 

GO  TO  262 
261  CONTINUE 

COMPUTE  THE  MAGNITUDE  OF  THE  SURFACE  FIELDS  AT  POINT  P  FROM 
CURRENTS  AT  POINT  PP. 

CALL  ESURFIPPX,PPY,PPZ,PX,PY,PZ, RADM, RADM, ETA, GAM , 

1  SPA , CPA , EMAGS , EMAGL . N ) 

IF  L  N  .  L  E  .  0  )  GO  TO  2S0 

COMPUTE  SHORT  PATH  AND  LONG  PATH  CONTRIBUTIONS  FROM  MODAL 
CURRENTS  1  AND  2  AT  LEFT  END  OF  THE  TEST  SEGMENT. 

DOT  CURRENT  DIRECTION  WITH  SURFACE  NORMAL  AT  SOURCE  POINT  PP. 

OOT-ABS(CAS*COSI  SPA)  )  +  ABS t CBS*S I Nl SPA  )  > 

Ell  1  “  EMAGS  +  EMAGL >*DOT*FAC*DK 

t.»2»(0.  ,0.  ) 


IF I NDK.NE .  1  )  En2«£NS 
262  CONTINUE 

CPA"ATAN2(PY,PX  ) 

DOT  THE  TOTAL  FIELD  FOR  THIS  OBSERVATION  POINT  WITH  THE  L-SEGMENT 
DIRECTION. 

ET1-EN1*<ABSICA*C0S(CPA>  )  ♦  ABS ( CB*S I N( CPA  !  > > 
ET2*EN2*(ABSICA*COSICPA>  )  ♦  ABSI CB*SI N( CPA  ))  > 

CC 

CC  ASSIGN  INTEGRATION  WEIGHTS  AND  MODAL  CURRENT  VALUES 
CC 

CW7-3 . +SGN 

IFIIN.EC.l  .OR.  IN.EO. IP >CWT»1 . 

EGD  =CE  XP I GAM*  I DL-T)  > 

C 1 *CVT*I EGD- 1 . /EGD>/2. 

EGD=C EXPIG AM*T > 

C2*CWT*( EGO- 1 . /EGD  > / 2 . 

P  1  1=P  1  1  +ET  1  *C  1 
P  12  =  P  12  +  ET1*C2 
P  2  1  =P2 i +E  T2*C l 
P22=P22+ET2*C2 
T=T+DELT 
100  SGN-  -SGN 

CST=  -DELT/I 3  .  *SG0T  ) 

PI  1 «CST*P 1  1 
P 12«CST*P  12 
P  2 1 =CST*P2 1 
P22*CST*P22 
P  I  1  .  I  )  ---  P  1  1 
PI  1  .  ?  >  =  P  1  2 
PI  2  .  1  >  =  P 2  1 
P i 2 , 2  )  =  P22 

168  C(I.n*C<I,J)+FI*Fvl*PIIS.OS) 

C  IF!  1  .  NE  .  J  >C!  J  .  I  >*C<  I  .  J  > 

200  CONTINUE 
250  RETURN 
E  ND 


Fig.  A-7.  (Continued) 


half  of  the  subroutine  is  the  same  as  for  SCANT G  or  SG.  However, 
unlike  SGANTC  or  SG,  the  current  on  the  source  segment  is  now  considered 
to  be  located  at  the  center  of  the  source  segment.  For  a  typical,  mid¬ 
antenna  segment,  the  current  distribution  is  made  up  of  two  modal 
currents  as  described  in  Chapter  IV.  The  magnitude  of  each  modal 
current,  previously  determined  in  the  main  calling  program  according  to 
the  procedure  detailed  in  Chapter  IV,  is  passed  to  subroutine  SGSURF 
through  variable  FAC.  Thus,  the  fields  ET1  and  ST2,  arriving  at  the 
various  observation  points,  P,  on  the  observation  segment  are  due  to 
current  moments  with  magnitudes  equal  to  FAC,  located  at  the  center  of 
the  source  segment,  generated  by  modal  currents  I  and  II  on  that  seg¬ 
ment.  Subroutine  VISIBLE  determines  if  the  source  segment  midpoint, 

PP,  can  be  "seen”  directly  by  the  observation  point  P,  The  diffracted 
field  contribution  is  determined  by  a  call  to  ESURF  only  if  the  two 
points  can  not  see  each  other,  ESURF  will  compute  the  values  for 
EMAGS  and  EMAGL,  EMAGS  is  the  field  which  arrives  at  the  observation 
point  due  to  the  shortest,  most  direct  path  between  the  source  and 
observation  segment,  as  shown  in  Fig,  4-3,  EMAGL  is  the  field  which 
arrives  at  the  observation  point  due  to  the  most  direct  path  between  the 
source  and  observation  segment  but  traveling  in  the  opposite  direction 
around  the  cylinder  as  Fig.  4-5  illustrates.  Subroutine  SGSURF  will 
then  add  the  two  fields,  EMAGS  and  EMAGL,  as  complex  numbers,  multiply 
the  result  by  the  length  of  the  segment,  and  take  the  necessary  dot 
products  as  required  by  eqn.  (3-1)  to  account  for  the  orientation  of  the 
source  and  observation  segments  with  respect  to  the  surface  of  the 
cylinder. 
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Subroutine  2SURF,  listed  in  Fig.  A-8,  implements  eqn,  (3-1). 

Based  on  the  location  of  the  source  point,  PP,  and  the  observation 
point,  P,  3SURF  determines  the  length  of  the  short  and  long  paths  around 
the  cylinder  by  applying  equs.  A-l  and  A-2  (see  Fig,  A-9).  The  pitch 
angle,  8  ,  is  derived,  based  on  eqn.  A -3  and  eqn.  A -4  provided  in 
Fig.  A-10.  The  subroutine  computes  the  other  parameters  of  eqn,  (3-1) 
in  a  straightforward  manner.  The  hard  and  soft  Fock  functions,  Fh 
(f,  ylt  7^)  and  Fs  (£,  y^,  7^)%  are  evaluated  by  calling  the  function 
subprograms  FH  and  FS  respectively. 

The  function  subprograms  FH  and  FS,  are  listed  in  Figs.  A-ll  and 

A-12,  To  increase  efficiency,  the  complex  values  listed  in  Table  3-1 

for  T  and  T  *  were  transformed  to  rectangular  notation  and  stored  in 
n  n 

data  arrays  TN  and  TPN.  The  constant  coefficients  for  the  hard  Fock 
function  residue  series  (eqns.  3-^  and  3-6)  are  stored  in  data  array  RK 
in  function  subprogram  FH  while  the  soft  Fock  function  residue  series 
(eqns.  3~5  and  3-7)  are  contained  in  data  array  RKK  in  function  sub¬ 
program  FS.  Data  array  KV  in  function  subprogram  FH  stores  the  constant 
coefficients  for  the  hard  Fock  function  small  argument  asymptotic 
expression  terms  (eqns.  3-9  and  3-H)  and  the  soft  Fock  function 
asymptotic  expression  coefficients  (eqns.  3-10  and  3-12)  axe  stored  in 
data  array  KU  in  function  subprogram  FS. 
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.0)8  RO’J  riNE  E  S  U  R  F  (  P  P  X  ,  F  P  v  .PPi.Pv.^V  ,  P  Z  ,  A  .  B  .  E  T A  .  GAM . SP A . C P A . 
1  EMAGS , EMAGL , N  ) 

COMP  .  E  X  ETA  ,  GAM  .  EMAGS  ,  LMAGL  .  EMAC- !  .  FMAG2 


C  *  NOTE:  A  A  8  >  N  THIS  SUBROUTINE 


ARE  IN  UNITS  Of  METERS.  NOT  * 
WAVELENGTHS.  * 

* 

C  «*»****»«»»*«*«««**»» 

c 

COMPLEX  FH.EHVAl ,  f S.ESVAL 
COMPLEX  J , G 
REAL  K 
REAL  M 

DATA  P I /3 . 14 1692653  3308/ . 0/1 «. . 1 .  )/ 

GiK\S;  =  K*CEXP(-J*'K''S>/<J*PI’*2.*S*ETA> 

IF ( A. NE . B  >  GO  TO  201 
SNMAG =SQRT(  *f‘F  X<  "PT  *P  PV  ) 

SPA  =  A  TA  N  2  (  I’PV  ,PPX  > 

C  A  MAG  =  SORT!  P  X  *  ’  X-*  P  V*P  V  > 

cfa-a r an r < pv  , p  » i 

D  1  --  SNMAG  -  A 

D2=CNMAG-A 

IF( SNMAG .LE .A  i  01=0. 


1 F ( CNHAG .1  E  .  A  )  DZ*0. 

C 

C  COMP  !  TC  1FNGTH  OF  SHORT  PATH  AROUND  THE  CYLINDER  FROM  NOPMAL 
C  PROJECTION  OF  POINT  P  TO  NORMAL  PROJECTION  OF  POINT  PP  ON  THE 
C  CYLINDER  SURFACE. 

i' 

S  >  =  SOP  T ( A*A*  C  S3A-CPA)**2.  ♦  <  PPZ-P  Z  )  *  *2  .  ) 

C 

C  SIMILAR. V,  COMPUTE  THE  LONGER  PATH  BV  GOING  THE  OTHER  WAV  AROUND 
C  THE  CYLINDER  BETWEEN  THE  TWO  POINTS. 

r_ 

SL=SORI(A*A*(  2.  “PI  -  AB3<  SPA -CPA  )  .'**2  .  ♦  (  P  P  Z- PZ  >  *  *  2  .  > 

< 

C  COMPurt  THE  PITCH  OF  TORSION  ANGi.ES  R OR  BOTH  SHORT  AND  LONG  PATHS. 

C 

li < PHZ-PZ.EQ.0. )  GO  TO  401 
PS IS=ATAN2< PP2-PZ,A*v  SPA-CPA  )  ) 

DEL  TS  =  ( PI/2.  )  -  APS! PS  IS  ) 

’SIL-/'TAN2(PPZ-P2.A*i<PI*2.  )  -  A6  S(  SPA-CPA  i  >  > 

>cl .  Tl  =  (  P  I  /.?  .  )  -  AES (  PS  I  L  ) 

GO  TO  40? 

401  l.‘ E  i.  T  !>  -  Pi/2. 

U  F  L  T  -  =  i  j  E  I.  T  S 
40.?  CONTINUE 
C 

C  COMPUTE  THE  MAGNITUCE  OF  THE  NORMAL  c I E  L  D  AT  P  FOR  BOTH  THE 
C  SHORT  AND  LONG  PATH  AROUND  THE  CYLINDER.  STARTING  WITH  THE  SHORT  PATH 

C 

>F  I.  T.w-E  L  TS 
C  -  s  s 

'  A  I  MAC ( GAM ) 

I'M  1.T0  [PATH-  I  .  2 
,  - 1  im  lA.tu.e  i  go  to  ?.pi 

V>  10  102 

.10  I  !  MAG  I  I  0.0.) 

Fig.  A-8.  Subroutine  ESURF 
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i  MAG  2  3  1  0  .  ,  0  .  ) 

GO  TO  303 

302  RHOG -A/(SlN(DtLTA>*SlNi DELTA)  ) 

RHOG *A/ ( SIN(DELTA)*SIN( DELTA)  ) 

M=( K “RHOG / 2 .  >**<  1  .  / 3  .  ) 

> I =M*S / RHOG 

!  F  (kI.LE.j3.)  GO  TO  904 
Y 1 -X  *0  1  /  M 
\  2-E*L2/M 
TO-COS ( DELTA  ) 
f HVA. =FH< XI  , Y1  ,  /2  > 

F  S VAL  =F  S ( X! , VI  .  Y2  > 

EMAG 1  - FH VAL  -  F  HVAL'O / <  K*S  >  +  f SVAL* ( 0 / ! K*S >  )* *  2 
EMAG  2  =  (  F  S  VAL  -  F  hVAI.  >  *  TO*TO*0  /  t  K*S  > 

303  IF  (I  PATH  EO.l)  GO  TO  201 
EMAG_=ETA*£TA*! t MAG  1  +  E HAG2 >*G ( X  ,  S ) 

GO  TO  100 

201  EMAGS=£TA*ETA*<  EMAG 1  +  E MAG 2  1  *G ( K  ,  S  1 
S=SL 

DELTA=OELTl 
100  CONTINUE 
RETURN 

ERROR  MESSAGES.  N  IS  USED  AS  AN  ABORT-COMMAND  FLAG. 

901  UR  I Tt  <  6 , 601  )  A , B 
GO  TO  310 

902  UR  ITS (6. 602)  SNMAG.A 
GO  M  310 


903  UP  IT:  (b, 603  >  CNMAG,  A 
GO  T  0  3  10 

904  UR  ITF  (  fj  .604  )  XI 
910  N - -  1 

RETURN 

601  FORMAT!//.*  Gi$$  FATAL  INPUT  DATA  ERROR*/ 

1  6*.  'GEOMETRY  EQUATIONS  IN  THIS  ROUTINE  ARE  VALID*/ 

2  6X , *ONL  V  FOR  CIRCULAR  CYLINDERS*/ 

3  SX.'A  -  * , E l T  3,ZX,*B  3  «,E10.3> 

602  FORMAT!//.*  SSSi  FATAL  INPUT  DATA  ERROR*/ 

1  6X  ,  *SOi.  RLE  POINT  IS  INSIDE  THE  CYLINDER*/ 

2  5X.*Si)MAG  *  *  ,E  10.3 . 2X  ,  *A  «  *.E10.3> 

603  FORMAT!//,*  Fatal  INPUT  DATA  ERROR.*/ 

1  6 < , 'OBSERVATION  POINT  IS  INSIDE  THE  CYLINDER*/ 

2  6a. ‘CNMAG  =  *.E’0.3.?X,*A  «  *.E10.3> 

604  FORMAT!//.*  $S$$  FATAL  INPUT  DATA  ERROR*/ 

1  b  X  ,  *  X  I  IS  LESS  THAN  OR  EQUAL  T0  0.0*/ 

2  6  X  ,  * X  I  =  * . E 10. 3  > 

END 

COMPLEX  FUNCTION  ES(XI,Y1,Y2> 

C  COMPUTES  THE  "SOFT'-TYPE  FOCK  FUNCTION 

COMPLEX  J  .  TN(  10  .  ,RKr-  <  2  > , XU ( 3  > . KUP (  3  > , SUM 
COMPLEX  ij,  up 
JA  TA  J  /  (  0  .  .  1  .  >  / 

DATA  R  K  K / ( 2 . 505628275  .2 . 30662327  5  >.(3.759942417,3.759942*12)/ 
DATA  TN/(  1 . 169055.-2 . 02 43631,(2. 043975, -3. 540269). 

:  (2.7  6028.-4.  7  809  4  5  )  . i 3.393  305 .  5. 877377), (  3.97206  5 .-6.879818), 

2  ( 4 . 5 1 1 325 , -7 .313844  ) . ( 5 .0201  . -8 . 59  5068 1,(5. 60425 ,-9.5336*1  )  , 

1  (5.965.-10  1316831  <6.4144.-11.1100671/ 

UAIA  K  U  /  (.  626')6  7068 /..  6/66570687  >.(  0  ...  4  1  66-666» />  , 

1  (.09/915167.  .097915167)7 

DATA  F  UP/< -.  939985003. -.  33998S603  .25). 
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I  ;  .  4406 1 025 t 4 , - . 4406 l 8 2 S 1 4  )  / 

IF  I  XI  .GT.J0.6  >  GO  TO  100 

COMF  ■  I.'  AM)  UP  USING  THt  FMALl  ARGlKtNT  ASYMPTOTIC  EXPRESSIONS. 

U-l.  -  KU(1)*aI“,*i.S  ►  K  j  t  *_  )  "  X  I  *  *  J  ♦ 

UP JP ( 1 )«SQRT. X 1 >  4.  K'JP  !  2  >*X  1  *  X  i  >  X  UP  '  3 X f *  *  3  S 

go  ro  ii0 

SPEAT’-P  THAN°ffU6  JSIH,S  ™E  C0NVERGING  RESIDUE  SCRIES  IT  XI  IS 

100  sum- ; 0. ,0.  > 

DO  200  N=  1  ,  10 

200  SUM-CEXPt -0»X1  *TN(  N  )  >  *■  SUM 
U i R  K K i  1  )  *  X  I  * ' 1  .  5  *  SUM 
S'JM  -  10.  ,0.  ) 

00  2  10  N  =  1  ,  10 

210  SUM*!  l"J«XI*TN<N>"2./3.  >  *  CE X P ( - J *X I  * TH ( N > >  -  SUM 
UP*P.*K.(  2  >*SORT<  XI  )  *  SUM 

C  COM* UTS  E  S ( X  I , V 1 , Y  2  )  AND  RETURN 

110  1  S  —  T  ►  J  *  .  5  *  (  i  J  P  -  1.5*"/Xr)*<y!*Yl+w2'Y'’> 

r-  I  ’.'IN 


(A-l) 

(A-2) 

(A-3) 


^  =  arctan 
s 


a  |  SPA  -  GPAJ 


|z2  “  zl| 

-  arctan  a'2^._  j  SPA  _  "cpA|) 

*s,./  =  f  '  *s,£ 


Fig.  A-lO,  Pitch  angle  calculations 
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I  ilMI'l  I  X  I  HNl  I  lull  I  M  V  1  .  V  |  .  v  .'  i 
l  I'l-H  HI  i  I  III  ••■,(i|  I  ”  I  VI'I  I  III  I.  IIINi  I  HIM 

I  '.I  in  i  X  .1  .  I  H'  1  I  >  .in  l.l  -•  I  ,  I  ll(  .1  )  .1  III’  (  !  )  .  SUM 
M.'l'i'M.l  X  U,  01* 

D  MA  J/(li.  ,  1  .  )/ 

IV,  f  \  KKK/<  2  .  ■  1  '.r.202  71,  .  r,0(,r>23'‘7‘l  )  ,  (  3 .  75904  .’41?  .  T  .  75994241  7  )/ 

L‘.\  1 IN/ (  I  1  ,.N''55 . -z  .0.'. I II 0  3  )  ,  <  2  .0439/ r. .  -X  510369  '  . 
i  <  2 . 7<>‘/?0  .  -  I  .  780945  ).(’■.  19 3 305 .  - 5 .  II  7  7  3  1  /  >  .  (  201.  •, .  T, .  17 'If, ;  R  ;  , 

:•  (  I  .  'i  1  I  '3  7  b  .-7.U13H44  >  .  I  5.0201  ,  -5  .  l.'lVMI  >  .  '  5  .  504  '5  ,  -•<).  5  ’Mil, 

3  ;  ,  .  .  -  li(.  3  3  I  5113  >.<(!..!  1  .HI  .  -  1  1  .  1  1 '  I  / 

DAI'/.  :.0/  (  .  0366570587  .  .  62665 7068 7  )  .  I  0 .  ,  .  -1 1  (.,,0  -6i.  17  >  . 

1  (  .9'.’  '315167. 097915167  )/ 

DATA  KUP/< 93998r-rj.il  1 ,  -  .  <1 3998560  3  )  .  ( IT.  ,  1  .  2  5  )  . 

1  (  .  I  4i,(il  020  l  4. -.4400  111251  4  >/ 

if  ( :;i  ,gi  .0.0 )  go  to  1.9.9 

CO;  IF  J  T  L  U  AND  UP  USING  THE  SMALL  ARGUMENT  AS  VHP  roil  t  r  '■ i  "F  '-Sinus  . 

Uu.  -  K  U  (  I  >■■*!{  I  *  *  1 . 0  +  KU< 2  )*XI**3  .  ♦  i:u<  3  Mill**  1  .0 
i"’-!;up(  1  >*suiu<  in  )  +  i;up<  2  >*xi*i:i  +  i:ur-<  3  >*::  1**3 .5 

GO  TO  1  19 

COMPUTE  U  AMD  UP  USING  1  HE  CONVERGING  RESIDUE  SERIES  IF  XI  IS 
GREAIER  TUAN  0.0 

1U1)  M'lMl  1.  ,9.  ) 

DO  ;;o  i  h- l .  i/I 

200  SUM  C.  UP  (  -  J  *X  I  *  1  IK  N  >  >  +  SUM 

i'  i.i  i  >  'xi*  + 1 . o  isun 

SUM-1  .  .if.  ) 

DO  ’1  :  ||=1  ,  10 

21.9  SUIIM  1-,1  ••XI*M'N<H>*2./3.  >  *  CE  XP  (  -0*X  I  *TN<  N  >)  +  S'JM 
Li  I  ,.i .  11  .ill  i  ll"  .0  "I 

COMPUir  FS1XI.Y1.V2)  AND  RETURN 
110  F  S  -  U  *  J*.5MUP  -  l.b*-U/Xn*(Yl*Yl*Y2*Y2> 

RE  TURN 
END 


Fig.  A-ll.  Function  Subprogram  FS 


8.1 


o  o  o  r>  o  r>  o 


function  f  h  <  <  -i.v- 

C  COHP*jT£  »  THE  "  HAR  U  ”  -  T V P 1  i  OCX  M  '".  MJN 

COMPLEX  J.  TFN  I  10  )  .  KK  i  ?  »  .  k'V  <  !  .  ‘'M3!  S"M 
COMPi.  t  X  V,V1 

DA  T  A  TPN7(  . 509 39b , - . 892  MJ8  ,  {  M  674095  , -2 . 8  1  301  5  )  . 

1  ( 2 . 41005. -4 . 174329  >,<  2  -  3816 If, .  -  b.  337563  ).<  3. 68609,-6.384495  >  . 

2  <4. 2442 45, -7. 351248), .4  76-  1"  .  -3. 257942),  (5. 26385, -9. 117256), 


3  < 5.73755,-9.9377^8  l,  it  19 24 10. 725551 . / 

DATA  KV/(  .  31332  35343,  .3,  ,.,.86:4  3  )  ,  (0.  ,  .  1  16666566/  ) , 

1  (  .01  71351542  .  .0171351  S'*?  > 

DATA  kVl  /(  .  6266570687  .  5  . 6.,o  76687  )  ,  <  0.  ,  .  5823  33333  3  > . 

1  <  .  1370812338,  -.  1 3 708  :  i  Ml  .! 

DATA  07(0.  , 1  .  )/ 

DATA  RK/(1.  253314137,  -1,  ./Mji;  137),  (2. 506628  275,2. 506  628275  )  7 
IF  ( XI .GT.0.6  )  GO  TO  100 
C 

C  COMMUTE  V  AND  VI  USING  THE  M>ali  \PGUMENT  ASYMPTOTIC  EXPRESSIONS. 

V  =  1  .  -  KV(  1  )  *■  X  I  *  *  1  .5  ^  K  r.  2  .  *Al  **-3  .  ♦  K'/<  2  >  ►  s  i  «  *  4 . 5 

Vl  =  !.  +  kVK  1  .5  -  ,  v  1<  2  )*XI  »*3  .  -  KVI(31*X!*'*4.5 

GO  TO  110 

COMPUTE  V  AND  VI  USING  THE  CONVERGING  RESIDUE  SERIES  IF  XI  IS 
GREATER  THAN  0.6 

100  S JM=<  0.  ,0.  ) 

DO  200  N  =  1  ,  10 

200  SUM*CEXP( -0*XI*TPN( N >  )7T’N< N >  +  SUM 

V  =  RK; 1  )  “SORT ( X!  )  “SUM 
SUM=<0. .0.  ) 

DO  210  N=  1  ,  10 

210  S UM -CE XP i -J*X1*TPN(N  i  )  +  SUM 
V1=RK< 2  )*XI*'l  . 5* SUM 

COMPUTE  FH1X1.V1.V2)  AND  RETURN. 

110  F  H = V  -  0*VI*( Yl'Yi  ♦  Y  2  *  V  2  )  / <  4  .  *  X  l  > 

RETURN 
END 


Fig.  A-12.  Function  Subprogram  FH 
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